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ABSTRACT 


This  thesis  supplies  a  statistical  and  economic  tool  for 
analysis  of  the  failure  characteristics  of  one  typical  piece 
of  equipment  under  evaluation:  a  cam-driven  reciprocating  pump 
used  in  the  submarine's  distillation  system.  Comprehensive 
statistical  techniques  and  parametric  modeling  are  employed  to 
identify  and  quantify  pump  failure  characteristics.  Specific 
areas  of  attention  include:  the  derivation  of  an  optimal 
maximum  replacement  interval  based  on  costs,  an  evaluation  of 
the  mission  reliability  for  the  pump  as  a  function  of  pump 
age,  and  a  calculation  of  the  expected  times  between  failures. 
The  purpose  of  this  analysis  is  to  evaluate  current 
maintenance  practices  of  time-based  replacement  and  examine 
the  consequences  of  different  replacement  intervals  in  terms 
of  costs  and  mission  reliability.  Tradeoffs  exist  between 
cost  savings  and  system  reliability  that  must  be  fully 
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EXECUTIVE  SUMMARY 


Today's  Submarine  Force  is  faced  with  many  new 
challenges.  These  include  an  ever- expanding  role  in  Joint  and 
Maritime  operations,  despite  a  diminishing  operating  budget 
and  a  shrinking  force  level .  Submarine  maintenance  is  one 
key  area  that  must  be  carefully  evaluated  to  ensure  maximum 
operational  readiness  while  working  within  the  confines  of 
budgetary  constraints.  A  Maintenance  Effectiveness  Review, 
MER,  is  currently  in  progress  to  evaluate  the  applicability 
and  effectiveness  of  all  aspects  of  submarine  maintenance. 
Any  comprehensive  evaluation  of  maintenance  practices, 
however,  requires  a  detailed  analysis  of  equipment  performance 
in  terms  of  failure  characteristics  and  life  expectancy. 

This  thesis  supplies  a  statistical  and  economic  tool  for 
analysis  of  the  failure  characteristics  of  one  typical  piece 
of  equipment  under  evaluation:  a  cam-driven  reciprocating  pump 
used  in  the  submarine's  distillation  system.  The  purpose  of 
this  analysis  is  to  evaluate  current  maintenance  practices  of 
time-based  replacement  and  examine  the  consequences  of 
different  replacement  intervals  in  terms  of  costs  and  mission 
reliability.  Tradeoffs  exist  between  cost  savings  and  system 
reliability  that  must  be  fully  understood  prior  to  making  any 
policy  decisions. 

The  analysis  examined  the  failure  characteristics  of  a 
sample  set  of  61  pumps  installed  between  May  1987  and  December 
1993.  The  data  consisted  of  140  failures,  of  which  14 


required  pump  replacements,  and  2020  total  months  of  pump 
operation.  The  primary  source  of  the  data  is  the  Navy's  3-M 
system.  Pump  failures  are  classified  as  repairable  or  non- 
repairable  and  thus  require  replacement.  The  probability  that 
a  given  pump  failure  is  repairable  was  evaluated  to  be  0.9. 
Additionally,  a  new  pump  will  experience  on  average  10 
failures  cycles  prior  to  being  replaced. 

Trending  analysis  indicated  an  increasing  failure  rate 
with  pump  age;  a  sign  of  system  wearout.  The  increasing 
failure  rate,  likewise,  indicated  that  a  time-based 
replacement  maintenance  policy  may  be  warranted.  A 
Nonhomogeneous  Poisson  Process  (NHPP)  was  chosen  to  model  the 
pump's  failure  characteristics.  Maximum  likelihood  was  then 
used  to  estimate  the  NHPP's  parameters. 

The  use  of  a  stochastic  modelling  allowed  a  comprehensive 
evaluation  of  the  current  maintenance  policy  of  time -based 
replacement  of  the  pump  at  a  periodicity  of  36  months  for 
Trident  submarines  and  60  months  for  non-Trident  submarines. 
The  pump  is  the  same  for  both  platforms.  The  optimal 
replacement  interval  based  entirely  on  minimizing  the  long-run 
average  system  costs  was  determined  to  be  111  months  for  an 
average  lifecycle  cost  of  $2226  per  month.  This  periodicity, 
however,  resulted  in  undesirable  mission  reliability.  The 
resulting  probability  of  the  pump  completing  a  mission  of 
three  months  just  prior  to  replacement  was  only  0.23  with  1.47 
expected  pump  failures  during  the  mission.  A  replacement 
periodicity  of  36  months  resulted  in  a  mission  reliability  of 
0.76  with  0.27  expected  failures  at  an  average  lifecycle  cost 
of  $3674  per  month.  The  60  month  periodicity  had  a  mission 
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reliability  of  0.63  with  0.46  expected  failures  at  an  average 
cost  of  $2611  per  month. 

The  results  of  this  analysis  indicate  that  a  60  month 
replacement  interval  may  be  acceptable  based  upon  reliability 
reliability  requirements.  The  periodicity  of  36  months, 
however,  may  be  excessive  and  warrant  extension  to  60  months. 
This  thesis  does  not  attempt  to  quantify  any  minimum 
reliability  requirements.  It  is  evident,  though,  that  a 
replacement  schedule  based  strictly  on  economic  considerations 
is  unsatisfactory  from  a  reliability  standpoint.  This 
illustrates  on  of  the  most  powerful  uses  of  stochastic 
modelling:  the  ability  to  predict  and  evaluate  the 
consequences  of  maintenance  policy  decisions  on  system 
performance . 

The  scope  of  this  thesis  has  implications  beyond  the 
operation  of  this  one  pump.  This  case  study  of  the  cam-driven 
reciprocating  pump  illustrates  the  type  of  analytical 
techniques  necessary  to  perform  a  comprehensive  evaluation  of 
a  shipboard  system's  performance.  The  goal  is  to  provide 
decision-makers  with  the  in-depth  statistical  foundation  to 
make  sound  decisions  on  maintenance  policy;  decisions  that 
directly  affect  the  readiness  and  ability  of  the  U.S. 
Submarine  Force  to  assume  its  expanding  mission. 


IX 


INTRODUCTION 


A.  BACKGROUND 

The  end  of  the  Cold  War  and  growing  concerns  over  the 
national  debt  have  forced  the  Navy,  and  specifically  the 
Submarine  Force  to  reevaluate  its  current  operational  policy. 
One  of  the  areas  being  closely  evaluated  is  that  of  Submarine 
Maintenance.  RADM  (sel)  R.E.  Frick,  NAVSEA  Deputy  Commander 
(Submarines  Directorate) ,  recently  directed  an  evaluation  of 
all  submarine  maintenance  performed  at  Depot  (shipyard)  , 
Intermediate  Maintenance  Activity  (IMA),  and  Ship's  Force 
levels.  This  Maintenance  Effectiveness  Review  (MER)  is  a 
continuous  process  to  evaluate  the  maintenance  needs  of  the 
Submarine  Force  while  working  within  the  fiscal  constraints  of 
a  shrinking  budget.  NAVSEA  PMS390,  Submarine  Monitoring, 
Maintenance,  and  Support  Program  Office  (SMMSO)  is  the  lead 
activity  in  charge  of  coordinating  this  review. 

The  review  process  requires  an  in-depth  review  of  the 
applicability  and  effectiveness  of  all  maintenance  actions 
performed  by  the  submarine  force.  The  goal  is  to  eliminate 
redundant  and  ineffective  maintenance  while  maintaining 
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optimal  materiel  readiness,  within  a  limited  budget.  The 
concepts  of  Reliability-Centered  Maintenance  are  being  used  to 
evaluate  the  time-based  and  "fix  when  fail"  policies 
historically  used  by  the  Navy. [Ref.l]  The  key  to  a  successful 
shift  in  strategies  is  the  ability  to  monitor  equipment 
conditions  and  make  reliable  maintenance  decisions  based  on 
real-time  data  and  predictive  analysis. 

The  goal  of  this  thesis  is  to  address  the  analytical  needs 
of  the  Submarine  Force  necessary  for  a  thorough  assessment  of 
current  maintenance  policies. 


B.  CURRENT  ANALYSIS 

The  SMMSO  organization  consists  of  a  staff  of  approximately 
one  hundred  personnel  based  in  Washington,  D.C.,  as  well  as 
monitoring  teams  stationed  at  each  submarine  base.  The  focus 
of  SMMSO 's  work  is  submarine  equipment  performance  monitoring 
and  maintenance  recommendations.  On  a  periodic  basis,  SMMSO 
site  teams  perform  various  performance  tests  on  a  submarine's 
systems.  These  tests  involve  everything  from  the  vibration 
monitoring  of  a  pump  to  oil  analysis  of  shaft  lubricating  oil . 
The  results  of  these  tests,  as  well  as  any  significant 
information  on  past  problems,  are  analyzed  by  a  specific 
engineer  who  is  SMMSO 's  system  expert  on  that  equipment.  This 
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information  is  used  to  predict  future  equipment  performance, 
make  recommendations  on  corrective  actions,  and  make 
maintenance  deferment  and  equipment  replacement  decisions. 

In  the  past,  the  scope  of  the  statistical  analysis  was 
focused  on  short-term  predictions,  i.e.,  refit  to  refit,  which 
was  at  the  time,  adequate  to  formulate  maintenance  policy 
given  the  large  Defense  Budget  and  high  priority  on  submarine 
readiness.  These  policies  have  been  very  conservative  in 
nature  and  possibly  more  restrictive  than  necessary.  In  view 
of  the  shrinking  budget,  these  policies  should  be  reviewed  to 
ensure  the  optimal  use  of  availible  resources. 

C.  PROBLEM  STATEMENT 

The  specific  area  of  analysis  to  be  addressed  in  this 
thesis  is  the  performance,  from  the  standpoint  of  reliability, 
availability,  maintenance,  and  replacement  needs  of  a 
particular  cam-driven  reciprocating  pump.  The  Class 
Maintenance  Plan  (CMP)  of  this  pump  is  currently  being 
evaluated  for  effectiveness  and  applicability.  Mr  Richard 
Youngk  is  the  SMMSO  systems  engineer  conducting  a  performance 
analysis  of  the  cam-driven  reciprocating  pump.  This  thesis  is 
in  conjunction  with  his  efforts  and  builds  upon  his  research 
and  analysis.  [Ref .2] 
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The  selection  of  this  particular  pump  for  the  study  is  in 
itself  not  significant,  but  the  pump  is  typical  of  the  systems 
being  analyzed.  The  analysis  presented  in  this  thesis  focuses 
in  a  broader  context  on  the  specific  questions  being  asked 
more  broadly.  This  thesis  proposes  a  process  or  method  to 
systematically  analyze  the  pump  data  presently  available.  It 
is  believed  that  this  procedure  can  then  be  adapted  for  use  to 
analyze  other  systems . 

Specific  questions  to  be  addressed  in  the  analysis  are: 

1.  What  is  the  expected  number  of  pump  repairs  and 
replacements  required  over  a  certain  time  interval? 

What  is  the  associated  expected  repair  cost  and 
downtime? 

2.  Does  the  failure  rate  change  over  the  life  of  the  pump? 

3.  What  is  the  optimal  equipment  replacement  interval? 

4.  Can  a  useful  stochastic  model  of  pump  performance  be 
constructed? 

D .  FAILURE 

Before  proceeding  further  with  the  analysis,  it  is 
necessary  to  define  what  is  meant  by  the  teirm:  fsilux's . 
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A  failure  is  classified  as  an  event  or  inoperable  state, 
in  which  any  item  or  part  of  an  item  does  not,  or  would  not, 
perform  as  previously  specif ied . [Ref . 3]  This  paper  further 
classifies  failures  as  being  either  repairable  or  non- 
repairable.  A  repairable  failure  implies  that  the  system  can 
be  returned  to  operating  specification  with  the  replacement  of 
only  a  small  fraction  of  the  system's  parts  and  in  a  short 
period  of  time.  Non- repairable  implies  the  pump  must  be 
completely  replaced,  or  else  an  extensive  overhaul  is 
necessary  to  restore  pump  operation  within  allowable  limits. 
In  the  context  of  this  paper,  pump  replacements  are  limited  to 
replacements  resulting  from  a  non- repairable  failure  and  not 
replacements  based  upon  any  other  criteria. 


E.  FAILURE  MODELS 

Two  basic  models  often  used  in  evaluating  repairable 
systems  are  the  Homogeneous  Poisson  Process  (HPP) ,  and  the 
Nonhomogeneous  Poisson  Process  (NHPP)  .  Both  the  HPP  and  NHPP 
are  counting  processes  used  to  model  the  number  of  component 
or  equipment  failures  occurring  over  the  life  of  the 
equipment.  Figure  1.1  shows  the  life-cycle  of  a  piece  of 
equipment,  with  ti  representing  the  time  to  the  i'"'’  failure, 
and  Xi  representing  the  time  interval  between  the  (i-l)®*'  and 
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failures,  i.e.  Xi  =  ti-ti.i-  To  be  more  explicit,  x^i 
represents  the  time  from  the  (i-l)®*'  to  i'^’’  failure  for  item 
(here  pump)  k.  Replacement  of  pump  k  means  that  the 
repairable- failure  generating  Poisson  process  starts  again 
from  scratch. 


Now  let  N(t)  equal  the  number  of  failures  occurring  prior 
to  pump  age  t.  Then  N(t)  is  considered  a  counting  process, 
{N(t) ,  t^O},  if  the  following  conditions  hold: 

1.  N(t)  ^  0. 

.  2.  N(t)  is  an  integer. 

3.  If  s  <  t,  then  N(s)  ^  N(t)  . 

4.  For  s  <  t,  N(t)  -  N(s)  represents  the  number  of  events 
occurring  in  the  interval  (s,t). 
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1.  Homogenous  Poisson  Process 

The  Homogenous  Poisson  Process  is  a  special  counting 
process  which  is  commonly  used  to  describe  systems  with  a 
constant  failure  rate  X.  More  precisely,  the  counting 
process,  {N(t),  t^O},  is  said  to  be  a  Homogenous  Poisson 

Process  with  rate  X,  X  >  0  if: 

1.  N(t)=0. 

2.  The  process  has  independent  increments,  i.e.,  the 
number  of  events  occurring  in  disjoint  time  intervals 
are  independent . 

3.  The  number  of  events  in  any  interval  of  length  t=Tj-Ti 
is  Poisson  distributed  with  mean  Xt 


i.e.  for  n  =  0,1,2,..., 


Pr[N{T.)-N{x.)’  n]  -  e' 


T  ,<T^,  n-0, 1,  .  .  . 


Furthermore,  the  conditional  probability  that  the  system  will 
survive  to  time  ij  given  that  it  is  operating  at  time  Xi, 
denoted  by  R(Ti,Tj),  is 


R{x  X  )=e 
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2 .  Nonhomogeneous  Poisson  Process 

The  Nonhomogeneous  Poisson  Process,  also  called  a 
nonstationary  Poisson  process,  is  a  counting  process  similar 
to  the  HPP  except  that  the  rate  function  is  not  constant,  but 
a  function  of  system  age,  t,  denoted  by  X(t)  .  The  Xi's,  i.e. 
the  times  between  equipment  failure,  are  not  necessarily 
identically  distributed.  This  mathematical  structure  can 
represent  an  increasing  failure  rate  over  time  resulting  from 
0quipment  wear.  Pump  data  indicate  such  behavior. 

Specifically,  the  counting  process,  {N(t),  t^O},  is  said 
to  be  a  Nonhomogeneous  Poisson  Process  with  age-dependent 
repairable  failure  rate  X(t)  if: 

1.  N(t)  =  0. 

2.  N(t)  has  independent  increments. 

3.  The  number  of  events,  i.e.  pump  failures,  in  the 
interval  (Ti,Tj)  is  Poisson  distributed  with 
mean  m(Ti,  Xj)  , 


'j 


(1.3) 
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The  probability  of  n  failures  in  the  time  interval 


(Ti,Tj)  is 


r.  rw/  V  X  1  -('"(V'#)) 

PrlNi-cp-Ni-c)-  n]  ‘  e  ’  ^  - ;-i— 


The  reliability  of  the  system  at  any  time,  t,  depends  on 
the  age  at  which  the  most  recent  failure  occurred.  Suppose 
that  a  failure  actually  occurs  at  time  t  =  ij,  where  Xj  is  the 
age  at  which  the  failure  of  the  incumbent  system  (pump) 
occurs.  Then  the  probability  that  the  system  will  not  fail 
for  at  least  age  ij.  +  x  units  of  time  is 


^  -(wC-t.T+x)) 


(1.5) 
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II.  MODEL  DEVELOPMENT 


A .  DATA 

The  data  used  in  this  study  consist  of  a  sample  set  of 
sixty-one  pumps  and  140  failures  oven  the  obsenvation  peniod 
beginning  in  May  1987  and  ending  in  December  1993.  The  final 
observation  was  either  a  failure  and  replacement  prior  to 
December  1993,  or  the  pump  was  last  observed  still  in 
operation.  This  sample  set  does  not  include  all  pumps  in 
operation,  but  a  subset  for  which  there  exists  reliable  dates 
for  pump  installation.  Additionally,  the  pumps  were  carefully 
screened  to  ensure  that  all  pump  replacements  resulted  from 
non- repairable  failures.  The  installation  and  failure  times 
are  rounded  to  the  nearest  month  and  henceforward  all 
references  to  time  and  age  will  be  in  months.  The  time  in 
service  is  adjusted  to  remove  any  inactivation  period  of  two 
or  more  months  in  length  from  the  pump's  operating  age.  The 
pump  itself  has  no  runtime  meter.  Engineering  logs  monitor  a 
pump's  operation  and  would  provide  a  more  accurate  measure  of 
the  actual  operating  time  for  each  pump.  However,  it  is  not 
feasible  to  accumulate  and  analyze  this  information  in  a 
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reasonable  time  frame.  Failure  and  replacement  data  originate 
from  the  Maintenance  and  Material  Management  (3-M)  system 
which  reports  routine  maintenance  via  OPNAV  Form  4790/K, 
Casualty  Reports  (CASREP) ,  the  Submarine  Maintenance, 
Engineering  Planning  and  Procurement  (SUBMEPP)  system,  and 
input  directly  from  the  Fleet. 

Mr.  Richard  Youngk  conducted  a  Failure  Modes  and  Effects 
Analysis,  FMEA,  to  identify  the  critical  failure  modes 
characteristic  to  the  pump  and  characterize  the  required 
repair  actions.  Information  on  the  types  of  failures,  their 
frequency,  the  repairs  requirements,  and  the  trends  for  the 
different  failure  modes  is  essential  elements  in  any  analysis 
of  equipment  performance.  This  thesis  focuses  specifically  on 
the  occurance  of  failures  and  the  required  repair  action. 
Failures  are  not  distinguished,  but  treated  as  one  entity. 
The  lifecycle  cost  in  terms  of  material  and  labor  is  also 
examined.  Appendix  A  is  a  copy  of  the  actual  data  analyzed. 


B.  ASSUMPTIONS 

Several  assumptions  were  made  regarding  the  pump  failure 

data.  These  assumptions  are: 

1.  The  submarines  in  the  sample  set  have  relatively 
similar  operating  cycles,  i.e.,  pumps  on  different 
submarines  will  face  the  same  operating  conditions 
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and  undergo  roughly  the  same  number  of  hours  of 
operation  for  the  same  time  period. 

2.  Every  component  failure  causes  equipment  failure. 

3 .  Failures  are  immediately  evident . 

4 .  Pumps  are  only  repaired  or  replaced  at  failure  and  not 
in  anticipation  of  failure.  In  actuality  this  is  not 
always  the  case,  but  the  data  set  was  cleansed  to  ensure 
only  replacements  related  to  failures  are  counted. 

5 .  Consecutive  failures  on  an  individual  pump  are 
independent . 

6.  A  repair  returns  the  pump  to  full  operation,  but  does  not 
necessarily  restore  it  to  a  "good  as  new"  condition. 

7.  Equipment  repairs  consume  no  appreciable  time. 

8 .  A  replacement  constitutes  the  installation  of  a  new 
pump  or  a  complete  overhaul  of  the  current  pump. 


The  accuracy  of  these  assumption  and  the  accuracy  of  the  data 
will  be  discussed  further  in  Chapter  IV. 


C.  ANALYZING  THE  DATA  FOR  TIME  DEPENDENCY 

One  of  the  initial  steps  in  the  analysis  is  to  determine 
if  the  data  are  compatible  with  an  increasing  failure  rate 
over  time,  i.e.,  with  aging  or  wearout .  Evaluating  the  time 
intervals  between  failures  provides  some  insight  into  this 
matter.  For  a  constant  failure  rate  over  time,  the  intervals 
between  successive  failures  should  be  very  similar;  with  all 
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interval  times  following  an  exponential  distribution  with  the 
same  mean.  Care  must  be  taken,  however,  to  ensure  all  data 
has  been  taken  into  account.  Invalid  conclusions  can  be 
drawn,  for  example,  by  comparing  the  mean  time  to  first 
failure  with  the  mean  interval  time  between  the  first  and 
second  failure  if  all  pumps  have  not  failed  at  least  twice. 
Here  the  existing  life  of  the  pumps  with  only  one  failure  is 
not  accounted  for  in  the  mean  for  the  second  interval. 

A  better  approach  is  to  compare  the  mean  interval  time 
between  successive  failures  for  pumps  with  like  failure 
numbers.  Table  2.1  displays  the  mean  and  standard  deviation 
(STD)  for  the  interval  times  between  failures  for  pumps  with 
at  least  two  failures,  with  at  least  three  failures,  and  so 
on.  In  all  cases  the  overall  trend  is  a  decreasing  interval 
between  failures  over  time,  indicating  an  increasing  failure 
rate.  The  correlation  coefficient  between  successive  failure 
times  also  provides  evidence  of  an  increasing  failure  rate. 
Table  2.1  also  shows  the  correlation  between  the  interarrival 
times  of  successive  failures.  The  correlation  between  and 
X2,  the  interval  times  for  failures  one  and  two,  and  between 
X2  and  X3  are  not  significant.  The  negative  correlation 
between  X3  and  X^  and  between  X4  and  X5,  however,  indicates  a 
decreasing  interval  time. 
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TABLE  2.1:  COMPARISON  OF  INTERVAL  TIMES  BETWEEN  FAILURES 


FAILURES 

Xi 

X, 

X3 

X, 

X5 

Xe 

SAMPLE 

NUMBER 

2  (MEAN) 

12.81 

10.25 

36 

(STD) 

7 . 97 

8 . 03 

3 

11.79 

8 . 88 

24 

8 . 05 

6 . 54 

8.78 

4 

9 . 93 

10 . 07 

8.67 

5.13 

15 

6 . 82 

7 . 01 

7 . 90 

4.77 

5 

9.22 

9.11 

7.22 

6.78 

5.67 

9 

4 . 76 

6 .33 

4 .61 

5.29 

4.57 

6 

12 .00 

7.50 

8 . 00 

5.50 

2.50 

3 . 00 

2 

1 

3.5 

5 

3.5 

1 . 5 

1 

CORRELATION 

(X,.  X,.J 

.  029 

.  066 

-  .374 

-  .  113 

1 . 0 

- 

SAMPLE 

NUMBER 

36 

24 

15 

9 

2 

- 

A  more  rigorous  statistical  test  for  distinguishing  between 
a  constant  failure  rate,  i.e.,  a  HPP,  and  a  monotonic  trend  is 
the  Laplace  Test .  The  Laplace  Test  is  based  upon  the  HPP 
property  that  given  n  arrivals  in  time  (0,F) ,  the  unordered 
times  of  arrival,  denoted  by  T^,  Tj,  .  .  .  ,  T^,  are  the  ordered 
statistics  from  an  independent  uniform  random  variable  on 
(0,F)  . 
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Therefore,  the  test  statistic  for  the  Laplace  Test,  U, 


has  approximately  a  standard  normal  distribution.  Bates 
(1955)  showed  this  approximation  was  adequate  at  the  5% 
significance  level  for  n  ^  4. [Ref. 4] 

A  slight  modification  is  required  if  the  system  is  observed 
until  a  specific  number  of  failures  occur.  In  this  reference 
case  F  =  Tn  and  the  modification  is 


Cox  and  Lewis  (1966)  showed  that  the  Laplace  test  is 
optimal  for  the  NHPP  with  rate 


-oo  <  a,  p  <  00,  2^  0 


(2.3) 
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in  testing  the  hypothesis  of  a  constant  rate,  i.e.,  Ho:p=0, 
against  the  alternate  hypothesis  of  a  trend, 

The  test  can  be  further  modified  to  examine  a  series  of  k 
independent  systems  with  the  same  value  of  P .  The  new  test 
statistic  is,  [Ref. 5] 


u  = 


E  E 

>1  M  2  J.1 


1  * 

— E  «/■/ 
12  7?  . 


(2.4) 


This  procedure  was  applied  to  13  pumps  meeting  the  criteria 
of  four  or  more  failures.  Table  2.2  gives  the  sample  set  used 
for  the  test.  The  evaluation  of  two  pumps,  number  57  and  60, 
required  the  modification  of  Equation  2.2  since  the 
observation  period  ended  with  a  failure  for  both 
pumps . 

The  resulting  value  of  the  test  statistic,  U,  was 
U  =2.016  whicr.  resulted  in  a  p-value  of  .04  for  the  two-tailed 
hypothesis  test .  Thus,  the  null  hypothesis  can  be  rejected  in 
favor  of  a  changing  failure  rate  at  the  conventional  5o 
significance  level.  Further,  the  positive  value  of  the  test 
statistic  indicates  an  increasing  failure  rate  over  time. 
This  test  included  only  13  pumps  of  the  original  data  set  of 
61  pumps.  It  is  reasonable  to  assume,  however,  that  this 
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TABLE  2.2:  LAPLACE  TEST  FOR  DATA  TREND 


LAPLACE  TEST  FOR  A  TREND  IN 
H„:  3  =  0,  H,: 

THE  FAILURE  RATE 

3  #  0 

PUMP 

FINAL 

OBSERVATIO 

N 

(MONTH) 

STATUS 

C=RUNNING 

R=REPLACED 

NUMBER 

OF 

FAILURES 

SUM  OF 

FAILURE 

TIMES 

48 

26 

C 

4 

79 

49 

56 

C 

4 

124 

51 

50 

C 

4 

120 

52 

60 

C 

4 

127 

53 

24 

C 

5 

82 

54 

29 

C 

5 

89 

55 

37 

C 

5 

83 

56 

45 

C 

5 

118 

57 

44 

R 

4* 

67* 

58 

53 

C 

5 

187 

59 

62 

C 

5 

186 

60 

37 

R 

5* 

127* 

61 

53 

C 

6 

168 

TEST  STATISTIC  VALUE  :  U  =  2.016 

P-VALUE  =  .043 

NOTES:  *  THE  TOTAL  NUMBER  OF  FAILURES  IS  MODIFIED  FROM 

TO  rij.i  FOR  PUMPS  57  AND  6  0  SINCE  THE  FINAL 
OBSERVATION  WAS  A  REPLACEMENT. 

**  PUMP  50  IS  NOT  USED  SINCE  THE  ADJUSTED  NUMBER  OF 

FAILURES  IS  3. 
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result  is  representative  of  the  entire  sample  set.  The 
Laplace  Test  was  in  fact  performed  on  all  pumps  with  at  least 
one  failure.  This  resulted  in  a  value  U  =  1.912  and  P  =  .06. 
Although  possibly  not  as  accurate  as  for  the  case  with  n  ^  4, 
it  does  seem  to  indicate  the  overall  trend  is  an  increasing 
failure  rate. 


D.  CHOOSING  A  MODEL 

Given  that  the  rate  of  failure  occurance  seems  to  increase 

with  pump  age,  the  next  step  involves  determining  whether  a 

mathematical  model  can  be  constructed  to  accurately  simulate 

the  occurrences  of  failures .  The  nonstationary  trend  of  the 

data  indicated  that  the  times  between  successive  repairable, 

i.e.,  non-fatal,  failures  were  not  identically  distributed. 

A  particular  Non -homogeneous  Poisson  Process  (NHPP)  has  been 

selected  to  model  the  occurrences  of  failure  over  the  life  of 

the  pump.  This  approach  was  chosen  based  on  the  fact  that 

the  NHPP  assumption,  besides  being  mathematically 
tractable,  is  a  good  representation  of  many  systems  in  the 
real  world  because  it  is  a  consequence  of  the  "minimal 
repair"  assumption.  "Minimal  repair"  implies  that  the 
repair  involves  the  replacement  of  only  a  small  fraction 
of  a  system's  constituent  parts . [Ref . 6] 
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The  following  parametric  failure  rate,  X(t) ,  was  used  to 


model  the  increasing  trend  in  failure  over  time: 


A(/)=e 


-oo  <  a,  p  <  ^  ^  0 


(2.5) 


This  form  of  failure  rate  was  introduced  by  Cox  and  Lewis 
for  the  analysis  of  failure  data  for  aircraft  air-conditioning 
equipment  in  1966.  [Ref. 7]  Ascher  and  Feingold  (1966)  used 
this  model  to  analyze  the  performance  of  submarine  main 
propulsion  diesel  engines.  This  model  has  the  advantage  that 
since  equation  2.5  is  positive  for  all  value  of  a  and  (3,  no 
nonlinear  restrictions  are  necessary  for  the  estimators  of 
these  parameters . [Ref . 8] 

Unlike  the  Cox  and  Lewis  situation,  however,  each  failure 
may  result  in  either  a  repair  or  a  replacement.  Therefore 
each  failure  is  subject  to  a  probability  of  being  repairable, 
denoted  by  p(ti,j),  and  a  corresponding  probability  of  being 
fatal  and  requiring  pump  replacement,  denoted  by  q(ti_j)  with 
q(ti,j)  =  1  -  p(ti,j)  .  In  fact,  some  repairs  may  be  relatively 
easy  while  the  platform  is  on  a  mission;  others  have  greater 
operational  impact  in  that  they  may  require  mission 
termination  for  repair  to  be  made. 
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E.  MAXIMUM  LIKELIHOOD  ESTIMATION 

Ths  par'a.mstBi's  cx  and  (3  havs  bsen  sstiniatad  using  thG  tnsthod 
of  niaxitnuni  likslihood.  Ths  tnaxiiuuin  likalihood  ©stiniats, 
abbraviatGd  as  MLE,  for  [3  is  found  first,  thGn  a  is  found. 
A  numGrical  solution  is  nGCGSsary;  thor©  is  no  sittipl©  closGd 
form  solution  for  p. 

BGfor©  dGriving  ths  MLE,  it  is  nscsssary  to  dsfins  ths 
following : 

k  =  indsx  for  ths  nurabsr  of  pumps  in  ths  sampls  sst 

j  =  indsx  for  ths  individual  pumps,  j  =  1,2,  ..k 

i  =  indsx  for  ths  failurs  numbsr  for  pump  j,  i  =  1,2, ..n^ 

Uj  =  ths  numbsr  of  failurss  for  pump  j 

fj  =  tims  of  ths  last  obssrvation  for  pump  j 

*Nots :  tims  is  squivalsnt  to  pump  ags  in  this  analysis* 
ti,j  =  tims  of  failurs  i  for  pump  j 

p(ti,j)  =  probability  that  ths  failurs  occurring  at  tims 
ti,j  is  rspairabls. 

q(ti,j)  =  probability  that  ths  failurs  occurring  at  tims 
ti j  is  not  rspairabls  and  thsrsfors  rsquirss 
pump  rsplacsmsnt .  q(ti,j)=  l-p(ti,j). 

Fj  =  indicator  variabls  indicating  ths  status  of  ths  last 
obssrvation 

Fj  =  1  for  non -rspairabls  failurs 
Fj  =  0  for  pump  last  obssrvsd  opsrating 
A(tj^j)  =  rats  of  occurrsncs  of  pump  failurss  for 
pump  ags  ti,j. 

Ths  failurs  may  or  may  not  bs  rspairabls. 

A(ti j)  =  ths  intsgratsd  valus  of  ths  failurs  rats  from 

pump  installation  to  ags  ti.j,-  ths  sxpsctsd  numbsr 
of  failurss  (rspairabls)  up  to  ags  ti,j 
6  =  vsctor  of  paramstsrs,  0  =  (a,  |3) 

L(e)  =  Combinsd  Likslihood  Function  for  pump  failurss 
L^,(0)  =  Conditional  Likslihood  Function  for  pump  failures. 

NOTE;  If  the  last  observation  is  a  failure,  then  fj  =  tnj,j. 


20 


Appendix  B  contains  the  complete  derivation  of  the  maximum 
likelihood  estimators  for  a  and  3  and  the  following 
explanation  will  only  deal  with  the  final  results  of  that 
derivation . 

The  combined  likelihood  function,  L(0),  for  all  pump 
failures  reduces  to: 


i(e> 

>=l  M  7=1 


rin/'(',) 

>1  »=i 


(2.6) 


Substituting  the  A(t)  described  by  Equation  2.5  results  in  the 
following  likelihood  function: 


1(6  >  e 


jA  i-\ 


(2.7) 


The  conditional  distribution  of  the  observations  given  n^ 
events  for  pump  j  is  formed  by  dividing  the  combined 
likelihood  function  by  the  marginal  probability  of  pump  j 
having  n^  total  failures.  Taking  the  logarithm  of  this 
conditional  distribution  produces  the  conditional  log 
likelihood  function. 
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log  PE"/  -E">*‘»8(^^'^^-1) 

>1  M  >1  >1  >1 


(2.8) 


The  maxitnizing  value  for  [3  is  then  obtained  by  maxitnizing 
log  [LcO)].  This  can  be  accomplished  by  setting  the 
derivative  of  the  Conditional  Log  Likelihood  Function,  with 
respect  to  [3,  equal  to  zero  and  solving  for  the  root,  p. 


<illogL/p)] 


J.1  (.1  j-i  P  >1  (l-e 


(2.9) 


Since  no  closed  form  solution  for  3  exists,  a  numerical 
method  is  necessary  to  solve  for  the  root.  First,  d[log 
L,,(6)]/d3  is  plotted  over  a  range  of  3  for  a  visual 
approximation  of  the  equation's  root.  Figure  2.2.  This  "best 
guess"  is  then  used  with  Newton's  Method  to  find  a  more 
precise  value  for  the  root.  The  resulting  maximizing  value  of 
3  is  p  =  0.02258. 
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MAXIMUM  LIKILIHOOD  ESTIMATION  OF  BETA 


Figure  2.2:  Visual  Approximation  of  the  MLE 
for  3 


A  similar  process  is  used  to  find  the  value  of  a.  The 
conditional  likelihood  function,  however,  cannot  be  used  to 
solve  for  a.  as  it  was  for  3  since  the  conditioning  process 
removes  a  as  a  parameter.  Therefore,  L(0)  must  be  used  to 
find  the  maximizing  value  of  a.  Taking  the  logarithm  of  L(6), 
differentiating  with  respect  to  a,  and  setting  the  equation 
equal  to  zero  results  in  the  following  expression: 
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(2.10) 


mi  .  y  _e_ly 

do.  >1  p  >1 

Solving  the  equation  in  terms  of  a  results  in  the  following 
expression: 

o-\n-y^ - .  (2.11) 

.>1 

Inserting  the  maximizing  value  of  p  =  0.02258  into 

Equation  2.11  produces  the  resulting  maximizing  value  of  a 
is  o  =  -3.189.  Appendix  C  contains  the  MathCad  3.1  program 
of  Newton's  Method  and  the  calculations  used  to  determine  the 
parameter  values  of  p  and  a  . 

The  resulting  maximum  likelihood  estimator  for  the 
failure  rate  of  the  pumps  is : 

i  0.  (2.12) 
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The  confidence  intervals  for  &  and  p  are  obtained  by 
inverting  the  observed  information  matrix  to  form  an  estimate 
of  the  variance- covariance  matrix  for  &  and  p  .  The 
unconditioned  likelihood  function,  Equation  2.7,  is  used  for 
calculating  the  confidence  intervals  since  conditioning 
removes  a  as  a  parameter.  The  resulting  95%  Confidence 
Intervals  are  given  in  Table  2.3.  Note  that  zero  is  not 
contained  in  the  interval  for  3,  giving  further  indication  of 
an  increasing  failure  rate.  Appendix  C  contains  the  Mathcad 
3.1  program  used  for  the  confidence  interval  calculations. 


TABLE  2.3:  CONFIDENCE  INTERVAL  FOR  MLE  OF  a  AND  p 


95%  CONFIDENCE  INTERVAL  FOR  &  AND  p 

PARAMETER 

LOWER  LIMIT 

MLE 

UPPER  LIMIT 

ALPHA 

-3.50 

-3.189 

-2 . 87 

BETA 

0.012 

0.02258 

0.033 

F.  EVALUATION  OF  MODEL  FIT 

The  use  of  any  parametric  model  to  quantify  a  system's 
behavior  must  also  involve  a  goodness  of  fit  analysis  to 
determine  the  adequacy  of  the  model .  Two  methods  are  used  to 
assess  the  suitability  of  the  afore  mentioned  NHPP  as  a  model 
for  the  failure  process.  One  method  is  graphical  and  the 
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other  is  a  formal  statistical  procedure  using  Pearson's  Chi- 
Squared  Test . 

The  first  method  compares  the  an  empirical  failure  rate 
taken  over  equidistant  intervals  to  the  model's  failure  rate 
for  the  corresponding  times.  The  interval  length  selected  was 
eight  months.  During  each  interval  the  total  pump  exposure 
was  calculated  as  the  sum  of  the  number  of  pumps  operating 
during  each  month  of  the  interval.  Similarly,  the  total 
number  of  failures  for  each  interval  was  determined.  The 
interval  failure  rate  is  simply  a  point  estimate  formed  by 
dividing  total  failures  by  total  exposure.  The  interval 
midpoint  was  used  to  calculated  the  model's  failure  rate, 
X(t)=e“^^^  Table  2.4  contains  the  calculated  rates  and  Figure 
2.3  is  the  corresponding  graph. 

The  model  closely  approximates  the  interval  failure  rate 
up  to  the  last  two  intervals,  i.e.,  intervals  49-56  and  57-64 
months.  This  region  indicates  a  decreasing  failure  rate.  One 
should  note,  however,  that  the  total  exposure  for  these 
regions  is  relatively  small  compared  to  the  other  intervals. 
The  small  exposure,  therefore,  may  not  be  representative  of 
the  overall  trend.  Further,  one  must  recall  the  nature  of  the 
equipment  being  evaluated.  The  system  is  a  cam-driven 
reciprocating  pump;  a  mechanical  system  already  noted  to 
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TABLE  2.4:  COMPARISON  OF  INTERVAL  FAILURE  AND  THE 

MODEL  FAILURE  RATE 


INTERVAL 

(MONTHS) 

TOTAL 

EXPOSURE 

(MONTHS) 

OBSERVED 

FAILURES 

INTERVAL 

FAILURE 

RATE 

MODEL 

FAILURE 

RATE 

1-8 

488 

23 

0.04713 

0.04529 

9-16 

438 

30 

0 . 06849 

0 . 05426 

17-24 

372 

29 

0 . 07796 

0 . 06500 

25-32 

277 

23 

0.08303 

0.07787 

33-40 

188 

15 

0.07979 

0.09328 

41-48 

152 

14 

0.09211 

0 . 11175 

49-56 

76 

5 

0.06579 

0 . 13388 

57-64 

28 

1 

0 . 03571 

0.16038 

Figure  2.3:  Observed  Interval  Failure  Rate  vs .  the 
Model  Failure  Rate,  X(t)= 


27 


exhibit  indications  of  wearout.  Thus,  the  system  is  highly 
unlikely  to  experience  any  reliability  improvement  at  this 
point  in  life.  The  system  engineer  at  SMMSO  independently 
reached  a  similar  conclusion  while  doing  a  parallel  study, 
i.e.,  that  the  observed  decreasing  failure  rate  late  in  life 
was  a  statistical  fluxuation,  resulting  from  sample  size,  and 
not  truly  typical  of  system  performance.  Of  course  it  is 
possible  that  pumps  on  some  vessels  will  actually  improve  with 
age,  up  to  a  point.  This  feature  remains  open  for  further 
investigation. 

The  Pearson  Chi-Squared  Test  is  a  formal  statistical  test 
for  goodness  of  fit.  Here  the  expected  numbers  of  failures 
under  the  NHPP  assumption  are  compared  with  the  observed 
numbers  of  failures.  The  same  eight  month  intervals  are  used 
with  the  exception  of  the  last  two  which  are  combined  to 
provide  at  least  five  failures  per  interval.  Table  2.5 
contains  the  results.  The  chi-squared  value  is  -  9.267 
with  7-3=4  degrees  of  freedom.  This  produces  a  p- value  of 
0.05  <  p  <  0.1.  Such  a  p-value  is  not  ordinarily  considered 
to  indicate  significant  departure  from  the  basic  hypothesis, 
in  this  case  the  model . 

Note  the  large  contribution  of  the  last  interval .  The  same 
arguments  as  above  can  be  made  for  this  deviation.  Even  so. 
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one  cannot  reject  the  possibility  that  the  data  follows  the 
NHPP  model  at  a  5%  significance  level.  Therefore,  the  model's 
representation  of  the  data,  although  not  perfect,  appears 
acceptable.  Note  further  that  acceptance  of  the  model  is  a 
conservative  step  in  predicting  future  failure  characteristics 
for  the  pump. 


TABLE  2.5:  CHI -SQUARED  GOODNESS  OF  FIT  TEST 


INTERVAL 

(MONTHS) 

TOTAL 

EXPOSURE 

(MONTHS) 

OBSERVED 

FAILURES 

EXPECTED 

FAILURES 

CHI 

SQUARED 

VALUE 

1-8 

488 

23 

22 . 1 

0.0341 

9-16 

438 

30 

23.8 

1.6172 

17-24 

372 

29 

24.2 

0.9468 

25-32 

277 

23 

21.6 

0 . 0910 

33-40 

188 

15 

17 . 6 

0.3735 

41-48 

152 

14 

17.0 

0.5324 

49-64 

104 

6 

15.3 

5.6717 

CHI  SQUARED  VALUE:  =  9.267  with  7-31  =  4  DF 

P-VALUE:  .05  <  p  <  .1 
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III.  DEVELOPING  A  REPLACEMENT  POLICY 


A.  BENEFITS  OF  STOCHASTIC  MODELING 

The  benefit  of  developing  a  stochastic  model  of  the  cam- 
driven  reciprocating  pump's  failure  characteristics  is  that 
pump  reliability  is  now  described  by  a  simple  mathematical 
expression  that  describes  likely  future  behavior.  One  might 
argue  that  empirical  data  would  provide  a  better 
representation  of  actual  pump  performance.  This  is  true  for 
systems  with  an  extensive  data  base  on  lifecycle  performance. 
Such  an  extensive  data  base,  however,  does  not  exist  for  the 
cam-driven  reciprocating  pump.  In  fact,  very  little 
information  exists  for  pumps  over  60  months  age.  The  lack  of 
such  data  demonstrates  the  need  for  a  probabilistic  model  to 
predict  future  behavior.  The  mathematical  model  allows  easy 
calculation  of  the  expected  number  of  pump  failures,  and  its 
availability  for  different  pump  ages  if  downtimes  are  also 
modeled.  It  also  expedites  determination  of  an  optimal 
replacement  interval  based  on  costs . 

This  chapter  will  discuss  the  merits  of  time-based 
replacement  and  evaluate  the  applicability  and  effectiveness 
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of  a  time-based  replacement  policy  for  the  cam-driven 
reciprocating  pump.  The  model  and  results  developed  in 
Chapter  II  provide  the  basis  for  the  evaluation. 


B.  CURRENT  MAINTENANCE  POLICY 

The  cam-driven  reciprocating  pump  is  currently  on  a  time- 
based  replacement  schedule.  The  schedule  varies,  however, 
between  platform  types,  specifically  between  Trident  and  non- 
Trident  submarines.  The  same  pump  is  used  in  both  platforms, 
although  slight  differences  in  the  systems  exist.  Current 
replacement  intervals  for  Trident  and  non-Trident  assets  are 
approximately  36  and  60  months,  respectively. 

The  time-based  replacement  policy  involves  the  replacement 
of  the  cam-driven  reciprocating  pump  at  a  predetermined  age 
regardless  of  the  current  material  condition  of  the  pump.  The 
value  of  this  time-based  replacement  in  maintenance  planning 
is  obvious.  A  planned  maintenance  evolution  allows  for  the 
scheduling  and  prepositioning  of  parts,  personnel,  and  support 
facilities  to  minimize  system  downtime  and  thus  tends  to 
minimize  total  system  (submarine)  nonavailability.  It  also 
precludes  possible  extensive  and  expensive  failures  that  may 
occur  at  later  ages . 
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This  reasoning  has  led  to  the  frequent  and  possibly 


excessive  use  of  time-based  maintenance  by  the  Navy . 

The  Navy  uses  time -based  maintenance  for  components, 
equipment  and  systems  ranging  in  complexity  from  oil 
filters  to  propulsion  gas  turbines.  Most  of  the 
maintenance  action  in  Class  Maintenance  Plans  are  based 
on  engineering  time  based  periodicities.  RCM 

[Reliability  Centered  Maintenance]  requires  that  these 
intervals  be  adjusted  based  on  equipment  performance  and 
failure  rate.  Most  time-based 

overhauls/refurbishments/replacements  are  also  expensive. 
CBM  [Condition  Based  Maintenance] ,  applied  where 
appropriate,  will  greatly  reduce  the  number  of  time-base 
repairs  and  overhauls  conducted.  The  key  to  successful 
implementation  of  CBM  is  application  of  the  proper  level 
of  monitoring,  evaluation  and  trending  for  each  piece  of 
equipment . [Ref . 9] 


The  challenge,  therefore,  is  to  evaluate  the  appropriateness 
and  effectiveness  of  the  current  time-based  maintenance 
practices  and  to  optimize  the  interval  used.  The  following 
issues  must  be  addressed  for  the  cam-driven  reciprocating  pump: 

.  1.  Is  this  pump  an  appropriate  candidate  for  time- 
based  replacement? 

2.  If  time-based  replacement  is  warranted,  what  is 
the  optimal  replacement  interval  to  minimize 
cost? 

C.  CRITERIA  FOR  TIME -BASED  REPLACEMENT 

The  goal  of  time -based  replacement  is  to  improve  both  the 
current  and  long-run  operating  state  of  the  system  through 
preplanned  maintenance  actions.  The  meaning  of  "operating 
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state"  is  dependent  upon  the  nature  and  mission  of  the  system  in 
question  and  upon  the  objectives  of  the  policy  makers.  Thus, 
"operating  state"  refers  to  operating  cost,  system  reliability, 
maintainability,  and  so  on. 

MIL-STD-2173 (AS)  provides  the  following  guidance  on  the 
applicability  of  time-based  tasks,  also  referred  to  as  hard  time 
tasks.  Here  the  time-based  task  is  considered  to  be  of  two 
types,  scheduled  rework  or  scheduled  discard.  A  reworking  task 
is  analyzed  if  reworking  promises  to  restore  the  item  to  an 
acceptable  level  of  failure  resistance;  otherwise,  the  discard 
task  is  analyzed. 

The  applicability  criteria  for  time-based  (hard  time)  tasks 
are  as  follows: 

1.  The  item  must  be  capable  of  having  an  acceptable 
level  of  failure  resistance  after  being  restored 
(for  rework  task) . 

2.  The  item  must  exhibit  wearout  characteristics,  which 
are  identified  by  an  increase  in  the 
conditional  probability  of  failure  with  increasing 
usage (age) .  This  property  can  lead  to  establishment 
of  a  wearout  age  (for  rework  tasks)  or  a  life-limit 
(for  discard  tasks) . 

3 .  A  large  percent  of  the  items  must  survive  to  the 
wearout  age  or  life-limit. 

4.  A  safe  life-limit  for  an  item  must  be  established 
at  an  age  below  which  relatively  few  failures  are 
expected  to  occur . [Ref . 10] 
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MIL- STD- 2 173 (AS)  points  out  two  key  elements  for  time-based 
replacement,  namely  an  increasing  failure  rate,  i.e.,  (2) 
above;  and  a  large  rate  of  survival  to  the  wearout  age  or 
life-limit,  (3)  above.  An  increase  in  the  occurrence  of 
repairable  failures  with  system  age  will  often  tend  to  result 
in  a  corresponding  increase  in  maintenance  costs .  A  decrease 
in  system  reliability  and  availability  will  also  occur.  From 
these  two  viewpoints,  optimal  replacement  intervals  should 
exist  to  minimize  costs  and/or  maintain  the  system  above 
certain  minimum  reliability  and  availability  requirements. 

A  large  percent  of  items  must  survive  to  the  point  of 
wearout,  i.e. ,  point  of  increasing  failure  rate,  or  life- 
limit,  to  make  planned  replacements  an  effective  maintenance 
tool.  Since  a  scheduled  replacement  is  deemed  to  be  more 
desirable  than  an  unscheduled  one,  the  opportunity  to  make 
replacements  should  be  utilized  fully.  Of  course,  the 
resources  associated  with  the  logistics  of  planning  a 
maintenance  evolution  may  be  wasted  if  a  large  percentage  of 
the  replacements  are  premature  and  occur  prior  to  the  planned 
interval . 

Several  additional  factors  must  also  be  considered  in 
determining  the  applicability  of  time-based  replacement.  Even 
a  system  with  a  constant  failure  rate  may  warrant  time-based 
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replacement  if  an  increase  in  operating  and  maintenance  costs 
occurs  as  the  item  ages,  or  if  repairs  at  sea  are  less  easily 
made.  Failures  may  occur  no  more  frequently  with  system  age, 
but  the  nature  and  type  of  repairable  failures  may  result  in 
an  increase  in  the  cost  of  parts  and  labor  as  the  system  ages . 
Another  consideration  in  a  time-based  replacement  scheme  is 
the  scope  of  work  required  to  replace  the  component  and  the 
maintenance  requirements  of  neighboring  systems.  Some 
replacements  may  require  extensive  interference  removal, 
elaborate  pre-established  plant  conditions,  and  extensive 
post- installation  testing.  In  such  cases,  common  sense 
dictates  combining  maintenance  actions  requiring  the  same  or 
similar  conditions.  Such  scheduling  may  not  coincide  with  the 
optimal  interval  to  minimize  the  operating  cost  for  every  item 
or  subsystem,  but  even  if  some  compromise  is  required  the 
overall  savings  could  be  substantial.  This  thesis  does  not 
address  the  problem  of  coincident  replacement  of  sets  of 
different  subsystems  that  have  age-dependent  failure 
properties . 


D.  APPLICABILITY  OF  THE  CAM-DRIVEN  RECIPROCATING  PUMP 

The  cam-driven  reciprocating  pump  is  a  candidate  for  time- 
based  replacement.  The  data  analysis  of  Chapter  II  suggests 
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an  increasing  occurrence  of  repairable  failures  as  the  pump 
ages,  thus  meeting  criteria  (b)  for  applicability.  As  stated 
earlier,  the  cam-driven  reciprocating  pump  is  currently  on  a 
time -based  replacement  schedule  of  36  and  60  months  for 
Trident  and  non-Trident  assets,  respectively.  The  remainder 
of  this  chapter  will  examine  available  data  for  evidence  that 
the  data  supports  the  current  replacement  intervals. 


E.  DERIVING  A  COST  MINIMIZING  FUNCTION 

As  Stated  earlier,  a  stochastic  model  of  the  pump's  failure 
characteristics  allows  the  use  of  mathematical  methodology  to 
predict  and  quantify  future  behavior.  One  such  tool  is  the 
Renewal  Reward  Process . 

Recall  that  the  Nonhomogeneous  Poisson  Process  is  used  to 
model  individual  pump  failure  times.  The  NHPP  is  not  a 
renewal  process,  but  since  all  new  pumps  are  assumed  to  be 
similar,  the  number  of  pump  replacements  to  occur  in  (o,t) 
does  constitute  a  renewal  process.  Let  M(t)  represent  the 
number  of  pump  replacements  occurring  in  a  system  up  to  and 
including  time  t,  and  let  L„,  n  ^  0,  represent  the  interval 
time  between  pump  replacements .  The  cost  associated  with  each 
renewal  is  denoted  by  R,,.  It  is  assumed  that  {L„}  and  {Rn}  are 
sequences  of  identically  distributed  random  variables;  a 
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generic  L„,  or  K,  is  denoted  by  L,  or  R.  Then  let  R(t)  be  the 
sum  of  all  system  costs  incurred  by  time  t,  so 

M(f) 

(3.1) 

n4> 

Further,  denote  the  expected  values  for  R„  and  L„  as  follows: 

E  [R]  =  E  [R„]  and  E  [L]  =  E  [L„]  .  Then  the  following  proposition 
holds : 

Proposition  3.1.  [Ref. 11] 

If  E  [R]  <  oo  and  E  [L]  <  «>,  then 


m] 

EIL] 


Let  a  cycle  denote  an  individual  pump's  life.  The  proposition 
states  that  the  long-run  average  cost  equals  the  average  cycle 
cost,  i.e.,  repair  and  replacement  costs,  divided  by  the 
average  cycle  length,  i.e.,  pump  life.  More  precisely,  the 
long-run  average  system  cost  equals 

E[cost  incurred  during  a  pump  's  life]  ^ 

1 

E[pump  life] 
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Pump  life  is  a  function  of  the  replacement  interval  and  the 
pump's  mortality.  Thus,  a  replacement  interval  for  minimizing 
costs  can  be  found  by  minimizing  the  long-run  average  cost,  as 
given  by  Equation  3.2,  over  different  replacement  intervals. 


1.  Expected  Pump  Life 

Let  X  be  a  random  variable  representing  the  age  of  the 
pump  at  replacement  and  let  T  be  the  designated  maximum 
replacement  age,  in  months.  Further,  define  a  cycle  to  be  the 
interval  from  pump  installation  to  pump  replacement  as  defined 
by  actual  replacement  or  complete  overhaul .  Then  the  cycle 
length,  denoted  by  L,  equals  X  if  the  pump  fails  and  requires 
replacement  prior  to  the  T;  otherwise,  L  equals  T  if  the  pump 
life  is  at  least  as  long  as  the  scheduled  replacement 
interval.  Thus  the  cycle  length  can  be  summarized  by 


X  if  OsX<r 
T  if  TiX 


(3.3) 


It  is  assumed  that  the  probability  a  failure  is 
repairable  is  constant,  i.e.,  independent  of  pump  age  and  the 
failure  number;  that  is,  a  failure  is  repairable  with 
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constant  probability  p.  The  data  have  140  failures,  of  which 
126  were  repairable.  This  results  in  a  point  estimate  for  p 
of  (126/140)  =  0.9.  The  approximate  normal  95%  Confidence 
Interval  for  p  is  (0.85,  0.95)  .  Likewise,  define  the 
probability  that  a  failure  is  not  repairable  as  q,  q=l-p. 
Brown  and  Proschan  used  a  similar  assumption  in  their 
imperfect -repair  model.  The  Brown  and  Proschan  model  assumes 
that  the  mode  of  repair  was  based  solely  on  external 
conditions  and  not  on  the  condition  of  the  system  at 
failure . [Ref . 12] 

Table  3.1  is  a  contingency  table  showing  the  distribution 
of  failures,  conditioned  on  failure  number.  The  assumption  is 
that  the  proportion  of  non- repairable  failures,  q,  remains 
constant  over  the  number  of  failures.  Let  rii  denote  the 
number  of  pumps  having  at  least  i  failures,  i  =  1,2,. ..,6. 
Further,  define  0,  as  the  number  of  observed  replacements 
and  Oi  l  as  the  number  of  repairs  for  pumps  with  at  least  i 
failures.  Define  the  total  number  of  replacements  and  repairs 
as  Ro,  and  R;,  respectively.  The  hypothesis  of  a  constant 
proportion  of  non- repairable  failures  is  tested  by  comparing 
the  observed  and  expected  values  under  the  null  hypothesis. 
Note  that  the  observations  for  failures  5  and  6  are  combined 
due  to  the  small  sample  size.  The  test  results  in  a  value  of 
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=  1.13  with  4  degrees  of  freedom  for  a  p-value  of  0.8  <  p 
<0.9.  The  hypothesis  of  a  constant  is  accepted  at  the  5% 
significance  level .  The  small  data  set  does  not  provide  for 
the  most  accurate  test,  but  it  does  give  some  indication  of 
the  goodness  of  fit.  So  the  assumption  of  a  constant 
probability  for  a  failure  being  repairable  is  not  unrealistic. 


TABLE  3.1: TEST  FOR  CONSTANT  REPLACEMENT  PROPORTION 


FAILURE  NUMBER 

STATUS 

B 

2 

3 

4 

mm 

TOTALS 

REPLACEMENTS 

5 

3 

2 

2 

2 

Ro  =  14 

REPAIRS 

49 

33 

22 

13 

9 

Ri  =  126 

bi 

54 

36 

24 

15 

11 

140 

CHI -SQUARED  VALURE :  =  1-13,  DF  =  4 

P-VALUE:  0.8  <  p  <  0.9 

M  ^  140  140j 

Now  define  h(x)  as  the  probability  density  function  for  a 
pump's  life,  X.  Then  h(x)  is 

K^)  -  E  ^  ”Hx)q  =  e  (3.4) 

f..o  «! 
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Similarly,  the  probability  that  the  pump  does  not  fail  before 
the  scheduled  replacement,  P(X  ^  T),  is 

p(x^Ty-  £  e -mWntp n  ^  (3.5) 

B-o  »! 

Using  the  above  information  it  is  possible  to  now  derive  an 
expression  for  the  expected  or  mean  life  of  a  pump  with  a 
designated  replacement  interval  of  time  T.  In  view  of  (3.3), 
the  expected  life,  denoted  by  E [L] ,  is 

T 

E[L]  -  jx  h(x)dx  .  fT  hix)dx.  (3.6) 

0  T 

Inserting  Equation  3.4  and  further  simplification  results 
in  the  following  expression  for  mean  pump  life: 

T 

E[Ly  jx  kix)qdx  .  T  (3.7) 

0 

The  above  general  formula  can  be  specialized  to  account  for 
any  parametric  form  for  A(t)  .  Often  the  formula  must  be 
numerically  evaluated;  closed-form  expressions  may  not  exist. 
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2 .  Expected  Cycle  Costs 

The  cost  incurred  over  the  life  of  a  pump  consists  of 
two  components;  the  cost  of  a  new  pump;  and  repair  costs 
incurred  until  the  next  replacement.  Thus  the  life  cycle  cost 
of  the  pump  is  influenced  by  the  number  of  repairable  failures 
occurring  over  the  pump's  life.  The  expected  life- cycle  cost 
can  be  represented  as 

E[life -cycle  cost]  =  COST  *  COST  ^^^E[number  of  repairs  ] 

where  both  COST„e„  and  COST^epair  are  expected  or  mean  values. 

As  stated  earlier,  the  pump's  life  can  terminate  in  one 
of  two  ways;  the  pump  may  have  a  non -repairable  failure  prior 
to  scheduled  replacement,  or  the  pump  will  be  replaced  per  the 
schedule  at  age  T.  Let  N(t)  represent  the  number  of 
repairable  failures  occurring  in  the  system  up  to  and 
including  time  t.  Likewise  let  E  [N]  represent  the  expected 
number  of  repairable  failures  occurring  during  pump  life. 
Then 


-  E[N{X)\X<T],P{X<n  +  ElNCD\XiT].P{XiT].  (3.8) 


The  analysis  will  examine  each  case  individually. 
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a.  Case  X  <  T 


The  expected  number  of  failures  given  that  the  pump 
does  not  live  to  the  scheduled  replacement  age  T  is 


5  [JV(r)|Ar<r]  ,P{X<T)  -  ^  nfe  K(x)qdx .  (3.9) 


This  simplifies  to 


T 

E[N(X)\}C<T\.P{X<T)  .  qpjA(x)e-^’=^k(x)dx.  (3.10) 

0 


Integration  by  parts  results  in  the  following  closed  form 
expression: 


£[Ar(r)lX<r].Pt!r<7’}  -  ^  -  ph(T)e-^  -  ^e'^. 


(3.11) 
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b.  Case  X  >  T 


The  expected  number  of  repairs  given  the  pump 
survives  to  scheduled  replacement  is 

E[N{T)\XiT\,P\XiT)  -  Y,  ■  (3.12) 

»0  "I 

This  expression  reduces  to 

E[N{T)\!C>^T].P{XiLT)  -  A(r)/»c‘*<^.  (3.13) 

Finally  combining  Equations  3.11  and  3.13  the 
expression  for  the  expected  number  of  repairs  over  the  pump's 
life  is 

^  -  pHT)e-^^  -  .  A(T)pe-^^.  (3.14) 

€  ‘I 

The  mathematical  expression  for  the  expected  cost 
incurred  during  the  pump's  life  is 

F[K].  COST  -  COST  -  pA(,T)e  -^^  -  .  A(T)pe-‘^'^ .  (3.15) 

"*"  npor  ^  J 


As  with  Equation  3.8,  the  above  general  formula  can  be 
specialized  to  account  for  any  parametric  form  for  A(t) . 


3.  Long-run  Average  System  Cost 

The  long-run  average  system  cost  can  now  be  expressed 
in  terms  of  Equations  3.7  and  3.15.  Let  z (T)  denote  the  long- 
run  cost  average  for  the  replacement  interval  of  length  T. 
Then 

-V-,  J. 

“  (3.16) 


Inserting  the  following  parametric  expressions  for  A(t) 
and  A(t) : 


A(x)  (eP*-l) 

P 

and  p  =  0.9,  the  optimal  replacement  interval  is  found  by 
minimizing  z (T) ,  displayed  in  Equation  3.16,  over  T. 


F.  DETERMINING  THE  OPTIMAL  PUMP  REPLACEMENT  INTERVAL 

The  first  step  is  to  determine  representative  values  for 
the  average  replacement  and  repair  costs .  Here  the 
replacement  cost  is  simply  the  cost  of  a  new  pump  which  is 


approximately  $111,000,  excluding  the  associated  installation 
costs.  The  repair  cost  is  a  function  of  failure  type.  The 
failure  type,  however,  may  be  affected  by  pump  age.  Some 
expensive- to- repair  failure  types  may  dominate  later  on  in  the 
pump's  life,  which  can  have  a  large  effect  on  the  average 
repair  costs.  Figure  3.1  exhibits  the  sample  distribution  of 
pump  failures  by  failure  type  as  the  pump  ages.  Note, 
however,  that  the  failure  data  is  not  adjusted  for  the  number 
of  operating  pumps.  Further  analysis  should  include  an 
estimation  of  failure  rate  for  the  specific  failure  types. 
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Figure  3.1:  Distribution  of  Failures  by  Type 
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This  leads  to  a  more  useful  representation  of  repairs 
cost:  a  moving  average  over  pump  age.  Using  information  on 
material  costs  for  repairs  from  the  3-M  System,  a  rough 
estimate  of  repair  costs  for  the  different  repair  types  was 
found.  The  3-iyi  data  is  sketchy,  however,  and  information  was 
not  available  for  all  repair  types . 

The  construction  of  the  moving  average  involved  ordering 
the  failure  times  and  assigning  and  an  average  repair  cost  to 
each  based  on  the  repair  type.  The  moving  average  consisted 
of  the  average  repair  costs  for  ten  consecutive  failures. 
Figure  3.2  show  this  average  cost  as  a  function  of  time.  Due 
to  the  limited  information  on  repair  costs,  no  attempt  was 
made  to  further  quantify  the  relationship  of  repair  costs  with 
age.  No  obvious  trend  is  evident  from  the  moving  average,  so 
for  a  conservative  estimate  of  repair  costs,  the  final  average 
occurring  at  month  61,  $4900,  will  be  used  as  a  basis  for 

replacement  interval  determination.  The  repair  cost  estimate 
deals  only  with  material  costs.  Other  factors  not  captured  in 
this  cost  estimate  include  labor  costs  and  overhead  costs. 
The  3-M  System  contains  limited  data  on  labor  manhours  for 
repairs,  but  again  this  information  is  sketchy.  Similarly  to 
the  material  cost  projections,  a  ten  point  moving  average 
was  constructed  to  assess  the  changing  maintenance 
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REPAIR  ACTION  MATERIAL  COST 


(10  FA.ILURE  MOVING  AVERAGE) 


Figxire  3.2:  Repair  Cost  Moving  Average  (10  pt .  Avg) 
requirements  as  the  pump  ages.  No  attempt  was  made  to  assign 
a  dollar  figure  to  these  labor  projections  due  to  the 
complexity  of  the  task  and  the  lack  of  quality  data.  Future 
research,  however,  should  be  directed  at  developing  such  a 
cost  relationship  to  be  combined  with  material  costs  for  a 
comprehensive  evaluation  of  trends  in  maintenance  costs. 

Figure  3.3  show  the  resulting  graph  of  the  moving  average 
for  expended  manhours  for  repairs  as  the  pump  ages .  The  graph 
shows  a  possible  increase  in  labor  requirements  with  pump  age. 
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As  stated  earlier,  no  attempt  is  made  to  incorporate  this  data 
into  the  cost  estimation. 
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Figure  3.3:  Moving  Average  Repair  Labor  Requirements 
(10  pt.  AVG) 

Determining  the  replacement  interval  is  now  simply  a 
matter  of  minimizing  z (T)  over  T.  A  Mathcad  3.1  program  was 
used  to  evaluate  and  graph  z  (T)  over  a  range  of  T.  The 
optimal  replacement  interval  is  the  minimum  point  on  the 
graph.  The  optimal  value  of  T  is  rounded  to  the  nearest 

month.  Table  3.2  contains  the  results  for  the  MLE  of  a  and 
3  as  well  as  for  their  95%  Confidence  Bounds.  Figure  3.4  is 
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a  graph  of  z (T)  for  the  MLE  of  a  and  p.  Calculations  are  made 
for  both  the  final  moving  average  value  of  $4900  and  the 
overall  average  repair  cost  value  of  $3000.  Calculations  are 
likewise  made  using  the  95%  Confidence  bound  derived  in 
Chapter  II. 


TABLE  3.2:  OPTIMAL  REPLACEMENT  INTERVAL 


C«e„  =  $111,000 

p  =  0.9 

REPLACEMENT  INTERVAL 
(MONTHS) 

Ot 

3 

c  = 

'-Repair 

$4900 

c  - 

'—Repair 

$3000 

-3.189 

.0226 

111 

12  8 

-3.50 

.  0226 

118 

135 

-2.87 

.0226 

104 

121 

-3 .189 

.  012 

183 

214 

-3.50 

.  012 

195 

227 

-2.87 

.  012 

171 

202 

-3.189 

.033 

82 

94 

-3.50 

.  033 

88 

99 

-2 . 87 

.  033 

77 

88 

The  results  indicate  an  optimal  replacement  interval  of 
111  and  128  months  for  the  average  repair  costs  of  $4900  and 
$3000,  respectively.  These  figures  as  well  as  the  calculated 
interval  using  the  confidence  bounds  are  well  above  the 
current  replacement  intervals  of  36  and  60  months  forTrident 
and  non-Trident  platforms,  respectively.  The  failure  and  cost 
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data,  therefore,  may  not  justify  the  current  Trident  and  non- 
Trident  platforms,  respectively.  The  failure  and  cost  data, 
therefore,  may  not  justify  the  current  replacement  interval  on 
the  basis  of  minimizing  maintenance  cost  alone. 


Figure  3.4:  Long-Run  Average  Cost 


Recall,  however,  that  only  material  costs  are  used  for  the 
repair  cost  estimate.  Other  costs  associated  with  repair 
include  the  cost  of  labor  and  shipyard  facilities.  Even  so, 
an  average  repair  cost  of  approximately  $32,000  would  be 
required  for  an  optimal  replacement  interval  of  60  months 
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using  the  MLE  parameters.  The  lowest  average  repair  cost 
required  for  a  60  month  interval  is  $11,000  and  occurs  at  the 
bounds  of  a  =  -2.87  and  |3  =  0.033.  Other  considerations  not 
accounted  for  here,  however,  may  be  included  in  the  current 
replacement  policy. 

In  developing  a  replacement  policy  for  the  cam-driven 
reciprocating  pump,  one  must  not  lose  sight  that  the  pump  is 
installed  in  a  warfighting  ship,  namely  a  submarine.  Cost  is 
not  and  should  not  be  the  lone  factor  in  establishing 
maintenance  policy.  System  reliability  and  its  effect  on 
overall  mission  accomplishment  must  be  a  strong  consideration 
in  any  decision.  To  aid  decision  makers,  the  following 
estimates  of  pump  performance  are  calculated  using  the  model : 
expected  number  of  failures  for  a  specific  pump  age,  the 
expected  failure  times,  and  pump  reliability  for  a  specified 


mission  duration. 

Since  a  NHPP  with  rate  X{t)  is  used  to  model  the  failure 
characteristics  of  the  pump  the  mean  value  function.  Equation 
1.3,  is  used  to  calculated  the  expected  number  of  failures 


over  pump  age . 


Figure  3 . 5  shows  the  expected  total  number  of 


failures  as  the  pump  ages  if  the  pump  is  never  replaced. 


Recall  that  each  failure  is  assumed  to  have  a  probability  of 
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q  =  0.10  that  the  failure  is  non- repairable .  The  sequence  of 
failures,  thus  constitute  a  geometric  distribution.  The  mean 
number  of  failures  before  a  required  replacement  is  (l/q)  or 
10  failures,  i.e.,  a  replacement  will  be  required  for  the 
tenth  failure. 


Figure  3.5;  Expected  Number  of  Pump  Failures 
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The  average  times  to  failure  can  be  calculated  using  a 
simple  simulation  program.  Since  the  NHPP  possesses  a 

continuous  mean  value  function  A(t),  (^inlar  (1975)  discussed 
the  following  recursive  algorithm  for  generating  a  sequence  of 
arrival  times,  t^tj/.-.t^: 

1.  Calculate  the  expectation  function  A(t) 
and  its  inverse  A“Mt)  : 

t  t 

A(t)  -  jk(x)dx  » 

0  0 

X  -  A(0  =  — (e  P'-l) 

P 

.  — hiC^+l) 

P 

2.  Generate  a  random  variable  U  ~  U(0,1) . 

3.  Set  t'i  =t'i.i  -  ln(U). 

4.  ti  =  A-^t'i]  .  [Ref  .13] 

Table  3.3  contains  the  mean  pump  age  at  failure  using  500 
replications  of  the  above  algorithm  for  simulated  failure 
times . 

The  decision  to  replace  or  repair  a  pump  may  be  guided  by 
certain  minimum  requirements  for  mission  reliability.  Recall 
from  Chapter  I  that  the  reliability  of  the  system  at  any  time 

t,  depends  on  the  age  at  which  the  most  recent  failure 

occurred.  Equation  1.5  gives  the  expression  for  the 

reliability  of  the  system  at  that  time.  Figure  3.6  is  a 

graph  of  the  probability  that  a  pump  having  a  repairable 
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TABLE  3.3:  SIMULATED  PUMP  FAILURE  TIMES 


FAILURE 

NUMBER 

MEAN 

PUMP  AGE 
(MONTHS) 

INTERVAL 

BETWEEN 

FAILURES 

STANDARD 

ERROR  OF 

MEAN 

1 

16.6 

- 

.602 

2 

29.2 

12.6 

.672 

3 

39.7 

10.5 

.691 

4 

48 

8.3 

.681 

5 

55.4 

7.4 

.659 

6 

61.8 

6.4 

.634 

7 

67.5 

5.7 

.624 

8 

72.4 

4.9 

.589 

9 

76.8 

4.4 

.555 

10 

81.1 

4.3 

.539 

11 

84.7 

3.6 

.517 

12 

88.1 

3.4 

.501 

13 

91.1 

3 

.489 

14 

94.1 

3 

.479 
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failure  at  age  t  will  successfully  complete  a  mission  of 
three  months  duration  without  experiencing  a  failure.  Figure 
3.7  shows  the  expected  number  of  failures  that  will  occur  over 
a  3  month  mission  given  the  pump  is  age  t  at  mission 
commencement.  This  is  computed  using  Equation  1.3  for  pump 
ages  ( t, t+3 ) . 

Both  Figures  3 . 6  and  3 . 7  indicate  that  the  replacement 
interval  of  111  months,  based  entirely  on  material  costs,  may 
not  be  desirable  due  to  the  poor  performance  of  the  pump  in 
terms  of  mission  survivability,  and  expected  pump  failures 
during  a  mission.  The  current  replacement  interval  of  60 
months  may  be  acceptable  based  upon  minimum  mission 
reliability  standards.  No  attempt  here  is  made  to  define 
those  standards.  The  replacement  interval  of  36  months  is 
probably  premature  and  should  be  extended  at  least  to  60 
months  to  coincide  with  the  policy  for  non-Trident  submarines. 

Table  3.4  provides  the  estimated  mission  survivability  and 
expected  number  of  failures  for  a  three  month  mission  for 
pumps  with  mission  completion  ages  of  36,  60,  and  111  months. 
Also  listed  is  the  long-run  average  cost  for  a  time-based 
replacement  schedule  with  the  associated  pump  ages .  This 
table  illustrates  the  trade-offs  in  cost  and  reliability  that 
must  be  resolved  for  an  effective  maintenance  policy.  All 
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Figure  3.6:  Reliability  for  a  90  Day  Mission 


Figure  3.7:  Expected  Number  of  Failures  During  a  90 


Day  Mission  as  a  Function  of  Pump  Age 
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decision  makers  should  understand  the  possible  consequences 


of  the  proposed  decisions.  The  use  of  probabilistic  modeling 
is  one  way  to  evaluate  the  trade-offs  and  consequences  of 
different  maintenance  policy. 


TABLE  3.4:  MISSION  SURVIVABILITY  AND  EXPECTED 


FAILURES  FOR  A  3  MONTH  MISSION 


PUMP  AGE 

AT  MISSION 

COMPLETION 

(MONTHS) 

PROBABILITY 

OF 

NO  PUMP 

FAILURES 

EXPECTED 

NUMBER 

OF  PUMP 

FAILURES 

LONG-RUN 

AVERAGE 

COST 

($ /MONTH) 

36 

.76 

.27 

3674 

60 

.63 

.46 

2611 

111 

.23 

1.47 

2226 

The  calculations  used  in  the  preceding  discussion  are 
relatively  simple  to  perform.  Many  similar  calculations  can 
also  be  made  to  address  specific  reliability  requirements  and 
different  measures  of  effectiveness.  The  goal  of  this 
analysis  has  been  to  provide  decision  makers  with  a 
comprehensive  assessment  the  cam-driven  reciprocating  pump's 
performance.  Such  information  is  necessary  in  the  formulation 
of  maintenance  policy  designed  to  address  both  economic  and 
reliability  concerns. 
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IV.  MODEL  ASSESSMENT 


A.  VALIDITY  OF  THE  DATA 

The  formulation  of  a  specific  NHPP  to  model  the  failure 
characteristic  of  the  cam-driven  reciprocating  pump  was 
accomplished  with  a  MLE  calculation  using  actual  failure  data,- 
data  obtained  primarily  through  the  Navy's  3-M  system.  The 
accuracy  of  the  model,  therefore,  depends  largely  upon  the 
degree  to  which  the  data  accurately  represents  actual  system 
performance . 

The  current  data  collection  method  is  far  from  perfect. 
The  3-M  system  itself  suffers  from  many  flaws.  Part  of  the 
problem  is  inherent  in  the  3-M  system,  itself,  and  part  is  due 
to  the  fleet's  attitude  toward  the  system.  When  a  failure 
occurs,  the  reporting  process  involve  a  crew  member,  normally 
junior  enlisted  personnel,  filling  out  an  OPNAV  Form  4790/K. 
Some  of  the  information  to  be  included  on  the  4790/K  are  the 
equipment  identification  code,  the  date  of  failure,  symptoms 
of  failure,  cause  of  failure,  required  repair  parts,  and 
required  repair  hours . 
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The  3-M  system  is  designed  to  track  failure  data  for 
everything  from  mechanical  and  electrical  systems  to  fuel 
storage  tanks.  In  doing  so  the  system  is  so  large  and  generic 
in  reporting  criteria  that  problems  arise  from  a  lack  of 
standardization  in  recording  specific  failure  information  for 
individual  pieces  of  equipment.  This  is  often  exemplified  in 
general  and  nonspecific  entries  for  failure  symptoms  and 
causes  which  can  lead  to  confusion  in  reconstructing  the 
actual  equipment  performance.  The  data  can  be  further 
confounded  by  improperly  entered  identification  codes  for 
equipment  and  repair  parts  as  well  as  incomplete  entries.  The 
3-M  system  also  does  not  capture  all  work  performed  by 
shipyard  personnel  during  non-availability  periods.  All  of 
these  factors  act  to  cloud  the  picture  of  true  system 
performance  and  thus  reduce  the  accuracy  of  any  data  analysis . 

The  data  used  in  the  formulation  of  this  model  may  not  be 
totally  accurate  in  its  portrayal  of  the  cam-driven 
reciprocating  pump's  maintenance  history,  but  it  is  currently 
the  only  viable  source  of  data  available  for  analysis. 

B.  VALIDITY  OF  THE  ASSUMPTIONS 

In  formulating  the  model,  eight  assumptions  are  made  as  to 
the  characteristics  of  the  cam-driven  reciprocating  pump 
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regarding  failures,  repairs  and  operation  of  the  pump.  Some 
of  these  assumptions  are  easy  to  accept;  others  require  some 
discussion.  First,  this  analysis  assumes  all  pumps  regardless 
of  the  submarine  in  which  they  are  installed,  experience 
roughly  the  same  operating  cycle.  Pump  operation  during  a 
submarine  deployment  will  be  similar  between  individual 
vessels.  The  deployment  schedule,  however,  will  differ 
between  submarine  platform  and  individual  units.  The  long-run 
average  deployment  time  is  assumed  to  be  the  approximately  the 
same  for  all  submarines  in  the  study.  This,  therefore,  is  a 
reasonable  assumption. 

The  assumptions  of  independence  between  consecutive 
failures  and  repairs  returning  equipment  to  full  operation  are 
related.  In  reality  these  assumption  are  not  always  true. 
One  failure  can  cause  a  subsequent  failure  at  a  later  age.  By 
the  same  token,  the  act  of  affixing  repairs  has  been  known  to 
cause  future  failures.  These  failures  may  be  totally 
unrelated  to  the  previous  failure,  but  were  caused  by 
improperly  restoring  the  equipment.  Likewise,  a  repair  may 
not  fully  repair  the  problem  and  the  pump  is  left  in  a 
condition  below  full  operating  capacity.  This  same  argument 
can  be  made  for  a  complete  overhaul  of  the  pump  in  that  the 
overhauled  pump  is  not  as  ' good-as-new' .  Postrepair  tests  and 
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procedures,  hopefully,  identify  and  correct  such  faulty 
repairs.  In  regards  to  this  data  set,  suspect  failures 
following  a  repair  were  carefully  scrutinized  to  catch  such 
double-failures.  In  the  end,  however,  all  such  dependence 
between  failures  and  incomplete  repairs  cannot  be  sifted  from 
the  data  set.  Therefore,  they  will  cause  some  loss  of 
accuracy  to  the  model . 


C.  ACCURACY  OF  THE  MODEL 

.beware  of  mathematicians  and  all  those  who  make 
empty  prophecies.  The  danger  already  exists  that  the 
mathematicians  have  made  a  covenant  with  the  devil  to 
darken  the  spirit  and  to  confine  man  in  the  bonds  of  Hell . 

St .  Augustine 

St.  Augustine  may  not  have  been  talking  about  modelling  the 
failure  characteristics  of  a  cam-driven  reciprocating  pump, 
but  he  does  provide  some  wisdom  for  modelling  in  general .  In 
fitting  a  mathematical  model  to  characterize  the  failure 
behavior  of  any  piece  of  equipment,  it  is  naive  to  think  that 
the  model  can  perfectly  predict  the  future  performance.  This 
would  be  true  regardless  of  the  system  or  the  quality  of  the 
data.  The  results  of  this  thesis,  therefore,  must  be  regarded 
in  this  light.  This  does  not  imply  that  any  such  model  has  no 
merit,  but  common  sense  and  good  engineering  principles  must 
be  incorporated  with  any  analytical  results .  The 
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incorporation  of  such  information  with  the  modelling  results, 
provides  a  sound  basis  for  making  policy  decisions. 

The  previous  sections  discussed  several  areas  that  may 
contribute  to  providing  inaccuracies  in  the  model,  namely  with 
the  accuracy  of  the  data  and  the  validity  of  the  assumptions. 
Another  factor  that  must  be  noted  is  that  the  original  data 
set  only  provides  information  on  pump  performance  up  to  about 
61  months.  The  majority  of  the  observations  actually  occurs 
below  40  months.  The  model,  however,  extrapolates  this 
information  into  performance  predictions  far  and  above  this 
age.  Any  model  inaccuracies  will  be  magnified  in  this  region. 
Also  the  use  of  this  particular  NHPP  is  only  one  parametric 
estimate  of  performance,  many  other  more  accurate  models  may 
exist. 

Given  all  the  model  inaccuracies,  does  this  model  provide 
any  information  to  policy  makers?  Yes!  The  formulation  of 
this  model  and  the  associated  data  analysis  have  revealed 
several  important  features  of  the  pump ' s  failure 
characteristics.  First,  the  pump  exhibits  wearout  as  noted  by 
an  increasing  failure  rate  with  age.  Secondly,  a  time-based 
replacement  decision  of  60  months  cannot  be  based  solely  on 
economic  reasons.  Thirdly,  the  model  provides  an  estimate  of 
the  expected  failures,  failure  times,  and  mission  reliability 
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as  the  pump  ages .  This 
precise,  is  pertinent 
maintenance  policy. 


information,  while  not  being  totally 
to  any  decision  maker  evaluating 
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V.  CONCLUSION  AND  RECOMMENDATIONS 


As  stated  in  the  Introduction,  the  goal  of  this  thesis  has 
been  to  address  the  analytical  needs  of  the  Submarine  Force 
necessary  for  a  thorough  assessment  of  current  maintenance 
policies.  Towards  this  end,  an  analysis  of  the  failure 
characteristics  of  a  cam-driven  reciprocating  pump  has  been 
conducted  to  demonstrate  the  degree  of  analysis  required  for 
a  comprehensive  assessment  of  equipment  performance.  Such 
analysis  is  necessary  to  assist  decision  makers  in  formulating 
maintenance  policy  especially  when  the  Navy  is  faced  with  the 
reality  of  a  shrinking  budget.  Decisions,  however,  must  not 
be  made  strictly  based  on  monetary  measures.  A  thorough 
understanding  of  the  consequences  of  any  decision  in  terms  of 
system  performance  and  reliability  is  essential. 

This  analysis  has  been  an  extension  of  the  work  preformed 
by  Mr.  Richard  Youngk  at  SMMSO  in  evaluating  the  operation  of 
the  cam-driven  reciprocating  pump.  His  assistance  and  the 
support  of  SMMSO  was  vital  in  the  completion  of  this  work. 

The  results  of  this  thesis  are  not  perfect.  This  analysis 
used  one  particular  model,  namely  a  NHPP  with  a  specific  rate 
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function,  to  model  failure  characteristics.  Other 
probabilistic  models  exist  and  a  different  model  could 
possibly  have  provided  a  more  accurate  representation  of  pump 
performance.  It  could  be  advantageous  to  compare  the  results 
of  this  thesis  with  some  of  the  other  models.  Additionally, 
one  slight  modification  to  the  current  model  is  recommended; 
perform  a  logarithmic  transform  on  the  failure  times  and 
recalculated  the  MLE.  This  would  provide  a  smaller  increase 
in  the  failure  rate  over  time  and  may  improve  the  model's  fit. 

As  stated  earlier,  the  accuracy  of  this  model,  and  in 
general  any  model,  is  to  a  large  degree  reflected  by  the 
quality  and  quantity  of  data.  The  modelling  process  is  not 
complete;  as  more  information  becomes  available,  the  model 
should  be  updated.  Likewise,  future  work  should  be  devoted  to 
improving  the  data  collection  system.  If  the  3-M  system  is  to 
provide  quality  information,  improvements  must  occur, 
specifically  in  the  timeliness  and  accuracy  of  reporting 
p0]ff ormance  data .  Many  independent  data  tracking  programs 
have  grown  out  of  frustration  with  the  current  system.  Such 
systems  while  providing  quality  information,  place  additional 
burdens  on  analyst  and  fleet  personnel.  The  3-M  system  needs 
to  re-examine  its  purpose  and  assess  how  well  it  addresses  the 
analytical  needs  of  the  Navy . 
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The  evaluation  of  the  cam-driven  reciprocating  pump  should 
continue.  New,  quality  data  will  improve  the  analysis  and 
better  quantify  pump  performance.  Specifically,  comprehensive 
cost  and  labor  estimates  will  more  clearly  define  the  changing 
maintenance  requirements  of  the  pump  as  it  ages . 

This  analysis  has  looked  at  only  one  small  aspect  of  pump 
performance.  Many  additional  areas  merit  research.  One  area 
centers  around  identifying  differences  in  pump  performance 
between  platforms  and  individual  submarines .  The  recognition 
and  investigation  of  such  differences  could  identify  specific 
operating  and  maintenance  practices  unique  to  certain 
submarine  that  either  enhance  or  degrade  system  performance. 

Another  area  for  further  study  involves  analysis  of  the 
failure  patterns  by  specific  failure  type.  The  thesis  model 
considered  failures  as  one  entity  when  in  fact  many  different 
failure  types  exist.  If  a  pattern  follows  the  occurrence  of 
certain  failures,  then  future  failures  and  failure  types  could 
possibly  be  predicted  based  on  previous  failure.  This  would 
aid  in  planning  maintenance,  but  more  importantly  it  would 
steer  investigation  toward  why  such  tendencies  exist  which 
could  lead  to  equipment  or  procedural  modifications. 


67 


The  concept  of  failure  prediction  and  prevention  through 
monitoring  also  deserves  attention.  This  again  involves  an 
analysis  of  failures  by  failure  type,  but  additionally 
reguires  the  identification  of  any  precursors  or  indicators  of 
impending  failure.  Such  monitoring  and  predictive  analysis 
already  exist  for  many  system.  Further  research  could  explore 
the  value  of  certain  precursors  including  determination  of  the 
expected  time  to  failure  given  a  specific  precursors  exist. 
This  could  aid  in  evaluating  the  effectiveness  of  monitoring 
procedures  and  possibly  provide  an  optimal  interval  for 
planned  monitoring. 

One  final  area  of  research  is  that  of  time-based 
replacements  for  multiple  systems.  This  concept  was  discussed 
briefly  in  Chapter  III  in  which  maintenance  for  different 
systems  may  be  scheduled  to  coincide  if  significant 
interference  removal  or  abnormal  plant  conditions  are  required 
for  both.  This  thesis  only  addressed  one  system,  but  multiple 
systems  could  be  analyzed  similarly  to  determine  the  most 
effective  replacement  schedule  to  maintain  the  readiness  of 
all  systems  involved. 

This  thesis  has  demonstrated  the  benefits  of  quantitative 
analysis  and  stochastic  modelling  in  evaluating  equipment 
performance.  Such  analysis  are  necessary  to  make  the  tough 
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decisions  on  maintenance  policy  in  the  wake  of  a  shrinking 
budget.  SMMSO,  and  specifically  Mr.  Richard  Youngk,  has  taken 
on  the  task  of  incorporating  these  and  similar  concepts  into 
current  maintenance  planning.  It  is  hoped  that  policy  makers 
will  use  these  mathematical  tool  to  make  the  well  informed 
decisions  necessary  to  maintain  the  U.S.  Navy  in  top  materiel 
readiness . 
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APPENDIX  A.  PDMP  FAILURE  DATA 


FINAL  EXCESS  TOTAL 


PUMP 

Ti 

T. 

T3 

T, 

Ts 

Tt  STATUS 

LIFE 

LIFE 

1 

c 

43 

43 

2 

c 

41 

41 

3 

c 

21 

21 

4 

c 

10 

10 

5 

c 

32 

32 

6 

c 

13 

13 

7 

c 

20 

20 

8 

1 

c 

9 

10 

9 

2 

c 

27 

29 

10 

2 

c 

22 

24 

11 

2 

c 

24 

26 

12 

3 

c 

11 

14 

13 

4 

c 

9 

13 

14 

5 

c 

22 

27 

15 

8 

c 

2 

10 

16 

10 

R 

0 

10 

17 

11 

R 

0 

11 

18 

13 

C 

1 

14 

19 

13 

R 

0 

13 

20 

15 

C 

8 

23 

21 

17 

R 

0 

17 

22 

18 

C 

7 

25 

23 

26 

C 

7 

33 

24 

27 

R 

0 

27 

25 

29 

C 

7 

36 

26 

8 

10 

C 

11 

21 

27 

4 

14 

C 

8 

22 

28 

10 

14 

R 

0 

14 

29 

1 

18 

C 

34 

52 

30 

17 

21 

C 

19 

40 

31 

19 

24 

C 

13 

37 

32 

22 

28 

R 

0 

28 

33 

23 

31 

C 

22 

53 

70 


FINAL  EXCESS  TOTAL 


PUMP 

Ti 

T, 

T, 

Tb 

Tg  STATUS 

LIFE 

LIFE 

34 

14 

34 

C 

16 

50 

35 

16 

38 

R 

0 

38 

36 

18 

41 

C 

18 

59 

37 

26 

61 

C 

4 

65 

38 

11 

12 

15 

C 

11 

26 

39 

4 

6 

21 

C 

15 

36 

40 

5 

20 

23 

C 

30 

53 

41 

10 

17 

26 

R 

0 

26 

42 

17 

32 

33 

C 

3 

36 

43 

36 

37 

41 

C 

3 

44 

44 

16 

23 

43 

R 

0 

43 

45 

17 

25 

43 

C 

1 

44 

46 

18 

24 

56 

C 

7 

63 

47 

2 

8 

9 

10 

R 

0 

10 

48 

15 

17 

22 

25 

C 

1 

26 

49 

15 

32 

38 

39 

C 

17 

56 

50 

27 

32 

37 

42 

R 

0 

42 

51 

3 

27 

42 

48 

C 

2 

50 

52 

4 

19 

52 

52 

C 

8 

60 

53 

6 

15 

17 

20 

24 

C 

0 

24 

54 

9 

12 

16 

26 

26 

C 

3 

29 

55 

4 

6 

14 

25 

34 

C 

3 

37 

56 

5 

16 

18 

36 

43 

C 

2 

45 

57 

4 

10 

25 

28 

44 

R 

0 

44 

58 

19 

31 

42 

46 

49 

C 

4 

53 

59 

12 

36 

43 

44 

51 

C 

11 

62 

60 

11 

22 

25 

34 

35 

37 

R 

0 

37 

61 

13 

17 

30 

32 

36 

40 

C 

13 

53 

LSG£i^ju  t 

T,  = 

FAILURE 

TIME, 

IN  MONTHS, 

FOR 

FAILURE  NUMBER  i 

C  INDICATES  PUMP  LAST  OBSERVED  OPERATING 
R  INDICATES  PUMP  WAS  REPLACED  AT  LAST  OBSERVATION 

NOTES : 

ALL  TIMES  ARE  ROUNDED  TO  THE  NEAREST  MONTH 
TOTAL  FAILURES:  140,  TOTAL  REPLACEMENTS:  14 
TOTAL  EXPOSURE:  2020 
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APPENDIX  B.  DERIVATION  OF  THE  MAXIMUM  LIKELIHOOD  ESTIMATOR 


Let 

k  =  index  for  the  number  of  pumps  in  the  sample  set 

j  =  index  for  the  individual  pumps,  j  =  1,2,  ..k 

i  =  index  for  the  failure  number  for  pump  j,  i  = 

Uj  =  the  number  of  failures  for  pump  j 

fj  =  time  of  the  last  observation  for  pump  j 

*Note:  time  is  equivalent  to  pump  age  in  this  analysis* 

ti,j  =  time  of  failure  i  for  pump  j 

p(ti,j)  =  probability  that  the  failure  occurring  at  time 
t,  j  is  repairable. 

q(ti,j)  =  probability  that  the  failure  occurring  at  time 
t,,3  is  not  repairable  and  therefore  requires 
pump  replacement.  q(ti,j)=  l-p(ti,j). 

Fj  =  indicator  variable  indicating  status  of  the  last 
observation 

F,  =  1  for  non- repairable  failure 
F  =  0  for  pump  last  observed  operating 
X(ti j)  =  rate  of  occurrence  of  pump  failures  for  pump  age 
ti  ..  The  failure  may  or  may  not  be  repairable. 
2\.(ti  j)  =  the  integrated  value  of  the  failure  rate  from 

pump  installation  to  age  ti,j;  the  expected  number 
of  f ai 2 ures (repairable)  up  to  age  ti  j 
0  =  vector  of  parameters,  6  =  (a,  |3) 

L(e)  =  Combined  Likelihood  Function  for  pump  failures 
L^O)  =  Conditional  Likelihood  Function  for  pump  failures 

NOTE:  If  the  last  observation  is  a  non -repairable  failure  then 
fj  =  and  F  =  1. 
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Consider  the  likelihood  function  for  one  pump,  i.e.,  the 


1/0).  e 


)  fl 


(Bl) 


Note  that  in  (Bl)  the  dependency  of  A,  X  and  p  on  parameters 
0  is  suppressed. 

For  the  case  with  the  last  observation  being  a  non- 
repairable  failure,  fj  =  t„j,j  and  Fj  =  1  so  expression  (Bl) 
reduces  to  : 


1/0).  « 


e 


-EA(r,y>A(^y-l)I 


Equation  (B2)  simplifies  to: 


(B2) 


i/0>  e  J 


PitnJ>) 


(B3) 
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Combining  the  observations  from  all  pumps: 


HO)-  e  '' 


v'’V/ 


■  PC.} 


(B4) 


To  obtain  the  conditional  probability  density  function,  or 
pdf,  of  the  observations  given  nj,  j  =  l,2,...k;,  failures, 
divide  L ( 0 )  by 


This  results  in  the  conditional  likelihood  function: 


4(e> 


ft-; 

>■1  --1 _ 

II  Mf/' 


(B5) 
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SPECIAL  CASE:  MODEL  FOR  INCREASING  HAZARD 

In  the  remainder  of  this  thesis  we  consider  the  following 
specific  model,  developed  by  Cox  and  Lewis  (1966),  pp.  45-54: 


Then  the  combined  likelihood  function,  L(6),  becomes: 


m  -  e 


Similarly  the  conditional  likelihood  function,  L<,(6),  is 


/"“Ur'IIy 

>1 


\[{e 


Here  0  =  (a,  0)  and  any  parameters  used  to  specify  p(ti,j). 
Taking  the  logarithm  of  Lc(6)  results  in 


log  1^(6  >  »5Ilog(«/)  'log  Pt 


>1  u  >1 
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To  find  the  Maximum  Likelihood  Estimator  for  3,  take  the 
derivative  of  log  Lc(6) »  equation  (B8) ,  with  respect  to  3, 
set  it  equal  to  zero,  and  solve  for  3-  The  resulting 
derivative  is 


dQos  i/e)) 

dp 


k  "j  k  k 

E  E  *ij  ^tE  -E 

j.\  M  P  j-l  >1  (1-e”^ 


P/a 


(B9) 


Graphing  Equation  (B9)  over  a  range  of  values  for  3  and 
visually  locating  the  zero  intercept  produces  a  rough 
estimate  of  the  maximizing  value  for  3-  This  initial 
approximation  can  be  further  refined  by  a  using  a  numerical 
method,  such  as  Newton's  Method,  to  solve  for  the  root.  This 
reduces  the  mathematical  complexity  by  eliminating  the  need 
for  deriving  a  closed  form  solution  for  the  root  if  a  closed 
form  solution  does  indeed  exists. 

A  similar  method  is  used  to  find  the  MLE  of  a.  The 
Conditional  Likelihood  Function,  Lc(e),  however,  cannot  be 
used  to  solve  for  a  as  it  was  for  3  since  the  conditioning 
process  removes  ol  as  a  parameter.  Therefore,  L(6)  ,  Equation 
(B6)  must  be  used  to  find  the  maximizing  value  of  a.  Taking 
the  logarithm  of  L(0) ,  differentiating  with  respect  to  a, 
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and  setting  the  equation  equal  to  zero  results  in  the 
following  expression: 


d[log  1(6)]  ^  _  Q 

du  j.i  P  j.i 

Solving  the  equation  in  terms  of  a  results  in  the  following 
expression : 


a  -  h 


PE", 

>1 

/-I 


(Bll) 


Inserting  the  maximizing  value  of  3  into  Equation  (Bll) 
produces  the  resulting  maximizing  value  of  a. 

CONFIDENCE  BOUNDS  FOR  a  AND  3 

The  confidence  bounds  for  the  estimation  of  a  and  3  are 
derived  using  the  Fisher  information  matrix  to  obtain  the 
asymptotic  variances  and  covariances  of  the  MLE ' s .  The 
Fisher  information  matrix  is  the  composed  of  the  negative 
second  partial  derivatives  of  the  sample  log  likelihood. 
Inverting  the  Fisher  information  matrix  evaluated  at  the  MLE 
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for  a  and  (5  produces  an  estimate  of  the  variance/covariance 
matrix.  The  variance  can  then  be  used  to  derive  the 
confidence  bounds  for  the  MLE's. 

The  partial  derivatives  are  derived  using  the  log  of 
equation  (B6)  and  are  as  follows: 

a^iogL(e)^ 

a^iogL(e)  a^iogi(fl)  -  (B13) 

aoap  apa«  p^  ^  P 

ie  ^-^>-1) .  — E  (/>  \  — E  <jfe  \  (B14) 

ap^  p*  >1  p^  >1  P  >1 
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APPENDIX  C .  MATHCAD  3 . 1  PROGRAMS 
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★♦★★★★★★★★★★★★★★★★★★★★★★★**’*********************************’*^********************************************* 


MATHCAD  3.1  PROGRAM  TO  COMPUTE  MLE  OF  BETA 
********************************************************************************* *****************^******* 

j  =0..  60  k^l..  150 

n.  -READ(NU)  tt.  =READ(TTU)  f.  :=READ(TFU) 


fl(p) 


THE  DERIVATIVE  OF  THE  CONDITIONAL  LOG  LIKELIHOOD 
PLOTTED  OVER  VALUES  OF  BETA 


n2  last(y)-  1 
n2  =  3 

f'W'" 

y  =0.02258 

^n2 


0.023 

0.0225818402 
y=  0.0225823081 
0.0225823081 
0 
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★★*★★**★*★*★*★***♦*★*♦★★*★★★****★★**★*****★★*★★**★★★*****★**★★★*★********★♦******♦*****************♦ 


MATHCAD3.1  PROGRAM  TO  DERIVE  THE  95%  Cl  FOR  MLE'S  ALPHA  AND  BETA 
★♦**★*♦*★★*★*★***★**★★*♦★**★★★★*★**♦★****★*♦**♦****************************♦*****************•*■****** 

j  -0,.60 

n,  =READ(NU)  tt.  :  =  READ(TTU)  f.  :=READ(TFU) 


Second  partial  derivatives  of  the  unconditioned  MLF 


daa(a,b)  :  =  (e  -  l) 

b  A-i 

j 


dbb(a,b)  — 
b^ 


alpha  =  -  3,189 
beta  :  =  .02258 

daa(alpha,beta)  =-139.441  /-daa( alpha, beta)  -  dab( alpha, beta)  \  ‘ 

dbb(alpha,beta)  =-1.257*10^  (-dab(alpha,beta)  -dbb( alpha, beta)/ 

dab(alpha,beta)  =-3.56*10^ 

I  0.026  -7.335*10“'*  \ 

MAT  = 

\-7.335*10'*  2.873*10  ^/ 


95%  CONFIDENCE  INTERVAL  FOR  ALPHA 

LLA  =  alpha  -  1 . 96  (MAT^  ULA  =  alpha  +  1 . 96-  (MAT^ 

LLA  =-3.504  ULA  =-2.874 

95%  CONFIDENCE  INTERVAL  FOR  BETA 

LLB  =beta-  1.96  (MATj  j) ^  ULB  =beta+  1.96  (MATj  ^ 

LLB  =0.012  ULB  =0.033 
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MATHCAD  3.1  PROGRAM  TO  DETERMINE  THE  OPTIMAL  REPLACEMENT  INTERVAL 


Expected  life  of  the  pump  (ELife): 

T  =  timed  replacement  interval  in  months 
p  =  probability  that  a  faiiure  is  repairable 

q  =  probability  that  a  faiiure  is  not  repairabie  and  requires  pump  repiacement 

i  =  1..200  T.  :=i 


a  =-3,189  b:=  0.02258  p  =- 


q  =  1  -  p 


b 


clA(t)  :=e‘‘+‘’'‘ 


ttf  =  case  in  virhich  pump  is  repiace  before  scheduied  intervai  T 
ttT  =  scheduied  replacement  at  time  intervai  T 


ttf.  =  qxe' 

Jo 


•dA(x)  dx 


ttT.  :=T.  e 
1  1 


f(x)  :=q  x-e'^^’'^'‘'  dA(x) 


ELife.  =  ttf. -h  ttT, 
1  1  1 


Expected  reward/cost : 

CReplace  =111000 


CRepair  =4900  /  .A/rVq  p  -A(Ti).q  p  -A(Ti)-q\ 

ECost.  :=  CReplace  +  CRepaii'  A ^T.j-p-e  -i- -  -  A^T.j-p-e  - --e 
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Reward  /  Renewal  Process:  Long  Run  Average  Reward  (z): 


ECost. 
z  _  1 

'  ELife. 


min^  -  10000000 

Zp  = 1000000000000 

j  -1..300 

mm.:  =  until[(z.-z^^,),(zj] 


n2  =  last(mm) 


z,  =2.22645726-10^ 

n2-  1 

z  =2.22641167-10^ 

n2 

z.,  =2.2264619-10^ 

n2+  1 


2192 


n2  =  lll 


MATHCAD  3.1  PROGRAM  TO  SIMULATE  PUMP  FAILURE  TIMES 

Simulation  of  failure  times  using  method  of  Law  &  Kelton,  1982,  pp.  509-510. 

OR]GIN=0  n  =500  a  =-3.189  b  =.02258 


-1..14 


Uj  j  -md(r) 


i  -1  ..n 


tPoj  =0 


‘P.,j  "‘PCi-DJ-'K".,]) 


CALCULATION  OF  MEAN  TIME  TO  FAILURE  FOR  FAILURE  i 

'  ' 


^£ail7j 


avg^=(  0  16.668  29.288  39.841  48.154  55.474  61.885  67.59  72.539  76.899  81.285  84.829 
(88.276)  (91.262)  (94.209) 
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CALCULATION  OF  THE  STANDARD  ERROR  FOR  THE  SAMPLE  MEAN  TIME  TO  FAILURE  i 

j  :  =  l..n 


samplestd. 


^(6illj-avg,)^ 

•5 

^(fail2j-avg2)^ 

j  ! 

j 

n~  1  J 

samplestd^  ■  ~ 

n  -  1 

samplestd^  - 


J _ 

n  —  1 


samplestd^  - 


]^(fail4j-avg4) 

j 

n~  1 


^(feilSj^avgj)^ 

j 

.5 

^(fail6j-avgg)2 

j 

n-  1 

samplestd^  ■  ~ 

n  -  1 

samplestd^ 


■5  1 

i 

^(fail8j-avgg)^ 

j 

j 

n-  1 

samplestd^  .  — : 

n  -  1 

samplesldg  - 


^(fail9j-avgg) 

J _ 

n-  1 


y  ^£aillOj  -  avgj 


samplestd , 


^(feilllj-avgj, 


samplestd , 


^(feill2j-avgj^^ 

samplestd “  — -  ■ — — - 


samplestd. ^  ^ 


^(feill3j-avgj3) 

J  ...  _ 

n  —  1 


^(feill4j-avg,4) 


samplestd . 


stderror  -  -  *  samplestd 

.('N/n)  J 


stderror' =(  0  0.602  0.672  0.691  0.681  0.659  0.634  0.624  0.589  0.555  0.539  0.517  0.501 

(0.489)  (0.479) 
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